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ABSTRACT 


We  consider  a  model  of  an  integrated  voice  and  data  network.  This 
model,  which  lends  itself  to  analytic  and  algorithmic  solution,  is  formu¬ 
lated  as  a  convex  optimization  problem.  The  objective  function  comprises 
two  types  of  functions:  the  congestion  cost  functions  and  the  rate  limi¬ 
tation  functions.  The  congestion  cost  functions  act  to  limit  the  average 
traffic  entering  into  the  network  to  values  that  would  not  cause  network 
congestion.  The  rate  limitation  functions  ensure  that  all  conversations 
are  fairly  treated.  The  model  can  be  generalized  to  solve  problems  of 
networks  which  handle  n  types  of  traffic  that  have  different  levels  of 
delay  sensitivity,  where  n  ^  2. 

A  joint  flow  control  and  routing  algorithm  is  constructed  which 
uses  short  term  average  information  on  the  network  utilization  to  set 
the  voice  packet  lengths  and  data  input  rates,  and  to  determine  the 
routes  for  each  conversation.  The  voice  packet  lengths  and  data  input 
rates  are  set  in  such  a  way  as  to  achieve  an  optimal  tradeoff  between 
each  user's  satisfaction  and  the  cost  of  network  congestion.  Additional 
protocols  are  specified  for  dealing  with  such  issues  as  congestion  avoidance 
and  control,  and  for  implementing  flow  control  on  a  more  dynamic  basis 
than  the  quasi-static  joint  flow  control  and  routing  algorithm  can  handle. 
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CHAPTER  I 


INTRODUCTION 


1.1  Motivation  For  The  Integrated  Network 

The  development  of  the  integrated  voice  and  data  network  concept 
has  resulted  from  a  combination  of  several  factors  including  the  following 

(a)  the  low  utilization  of  network  resources  when  voice,  data, 
fascimile,  etc.  are  transmitted  on  separate  communications  networks  in 
which  traditional  switching  techniques  are  used; 

(b)  the  success  of  such  data  communications  networks  as  the  Telenet 
[33],  the  Tymnet  [30],  the  French  Transpac  [10],  and  the  Canadian  Datapac 
[7]; 

(c)  the  recent  technological  advances  in  speech  processing,  parti¬ 
cularly  speech  digitization. 

In  a  telephone  conversation  a  speaker  alternates  randomly  between  the 
speaking  (or  talkspurt)  mode  and  the  listening  (or  silence)  mode. 
Measurement  studies  have  shown  that  in  a  typical  telephone  conversation 
a  speaker  uses  the  channel  only  35-40%  of  the  time  duration  of  the  con¬ 
versation  [4].  This  then  means  that  dedicating  a  channel  to  a  pair 
of  speakers  throughout  the  duration  of  their  conversation  is  a  waste  of 
the  network  resources. 

Many  techniques  have  been  proposed  to  minimize  this  waste.  These 
techniques  attempt  to  compress  n  voice  conversations  onto  m  channels, 
where  n  >  m.  The  earliest  of  these  techniques  is  the  Bell  System  "Time 


Assignment  Speech  Interpolation"  (TASI)  in  which  a  channel  is  allocated 
to  a  subscriber  only  when  the  appropriate  hardware  detects  that  he  is 
actively  speaking  [4].  Once  he  accesses  the  channel  the  speaker  uses 
it  without  interruption  for  the  duration  of  his  talkspurt.  During 
periods  of  silence  he  relinquishes  the  channel,  and  it  becomes  available 
to  other  speakers.  Any  talkspurt  that  is  generated  when  all  channels 
are  busy  is  "frozen  out"  and  will  have  to  wait  on  a  first-come  first- 
served  basis  for  an  available  channel.  This  freeze-out  results  in  the 
clipping  and  loss  of  the  initial  part  of  the  talkspurt.  It  is  even 
possible  to  lose  an  entire  talkspurt  if  the  waiting  time  is  longer  than 
the  spurt  duration  [30].  For  a  satisfactory  speech  quality  the  fraction 
of  the  lost  speech  should  be  less  than  0.5%.  The  ratio  of  the  number 
of  subscribers  supported  by  the  TASI  system  to  the  number  of  channels 
required  to  maintain  an  acceptable  fraction  of  lost  speech  is  called 
the  "TASI  Advantage".  If  p  is  the  fraction  of  time  that  a  typical 
speaker  is  in  talkspurt,  the  maximum  potential  TASI  advantage  is  1/p. 
Since  p  <  0.5,  this  advantage  is  greater  than  2.  However,  at  least 
40  channels  and  80  speakers  are  required  to  achieve  this  TASI  advantage 
[4].  For  a  smaller  number  of  channels,  attempts  to  achieve  a  TASI  advan¬ 
tage  close  to  1/p  will  result  in  unacceptably  high  lost  fractions  [32]. 
The  TASI  system  described  above  is  oriented  toward  analog  speech  trans¬ 
mission.  Digital  variations  of  TASI  that  are  transmitted  in  the  PCM 
form  have  also  been  implemented.  These  include  Digital  Speech  Inter¬ 
polation  (DSI)  [5],  and  Speech  Predictive  Encoding  (SPEC)  [28]. 

Speech  is  bursty  in  nature,  and  this  accounts  for  the  low  utilizati 
of  the  transmission  resources  in  a  voice  communications  network  in  which 
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the  traditional  frequency  division  multiplex  (FDM)  or  time  division 
multiplex  (TDM)  is  used.  The  purpose  of  the  TASI  schemes  is  to  achieve 
a  better  utilization  of  the  network  resources  by  smoothing  out  the 
voice  flow.  Data  are  also  bursty  in  nature.  An  effort  to  smooth  out 
data  flow  in  a  network  has  led  to  the  concept  of  statistical  multiplexing 
[8]  or  store- and- forward  networks  [1].  Unlike  the  TDM  system  traditionally 
used  for  data  transmission,  no  slots  are  dedicated  to  any  data  source. 
Instead  messages  are  merged  into  a  buffer  at  each  network  node  and  the 
bits  are  transmitted  at  a  synchronous  rate.  The  messages  work  their 
way  through  the  network  to  their  destinations  .queueing  at  each  node  where 
they  are  statistically  multiplexed.  The  messages  may  remain  in  complete 
form  (message  switching)  or  be  split  into  shorter  blocks  called  packets 
(packet  switching)  [1],  Whether  statistical  multiplexing  exists  in  the 
message-switched  or  packet-switched  form,  each  unit  carries  a  header 
which  specifies  among  other  things  the  source  and  destination  addresses 
of  that  unit.  There  is  basically  little  difference  between  packet 
switching  and  message  switching.  We  shall  use  the  term  packet  switching 
hereafter  for  both  techniques. 

That  speech  is  bursty  in  nature  makes  it  a  good  candidate  for 
statistical  multiplexing.  This  has  been  further  made  possible  by  the 
recent  advances  in  speech  processing  which  have  led  to  the  packetization 
of  speech.  This  means  that  packets  of  speech,  which  are  generated  only 
when  a  speaker  is  in  talkspurt,  could  be  stored  temporarily  at  each 
node  in  wait  for  free  channels.  Thus  there  would  be  no  loss  of  any  part 
of  the  speech  as  is  often  the  case  in  the  conventional  TASI  system. 

And  the  fact  that  speech  and  data  can  be  packetized  raises  the  possibility 


of  using  the  same  network  to  transmit  them.  One  then  hopes  that  such 
an  integration  of  speech  and  data  onto  the  same  network  would  lead  to 
better  utilization  of  the  network  resources  [17,18]. 

1 .2  Problems  With  Integration 

The  integration  of  voice  and  data  traffics  onto  one  network 
generates  new  communications  problems.  These  problems  include  the  choice 
of  an  appropriate  switching  technique,  and  how  to  control  the  network 
flow  in  order  to  meet  the  different  demands  made  on  the  network  by  voice 
and  data  conversations. 

Three  switching  strategies  have  been  proposed  for  use  in  the 
integrated  network  [17].  These  are: 

(a)  Circuit  Switching 

(b)  Packet  Switching 

(c)  Hybrid  Switching. 

In  circuit  switching,  an  end-to-end  connection  is  established  for  a 
pair  of  voice  or  data  users  before  they  commence  their  conversation. 

The  channel  is  dedicated  to  them  throughout  the  duration  of  the  conversa¬ 
tion.  The  channel  utilization  may  be  improved  by  using  the  TASI  system. 
In  packet  switching,  both  voice  and  data  conversations  are  digitized 
and  segmented  into  packets  and  routed  through  the  network  in  a  store-and- 
forward  manner  [1,21].  That  is,  each  packet  travelling  from  source  node 
i  to  destination  node  j  is  "stored"  in  queue  at  an  intermediate 
node  k  while  awaiting  transmission.  It  is  sent  "forward"  to  a  selected 
adjacent  node  Z  on  the  route  to  node  j  when  link  (k,Z)  is  free.  The 
selection  of  node  l  is  made  by  some  well-defined  decision  called  the 
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routing  algorithm.  The  process  of  storing  and  forwarding  is  repeated 
until  the  packet  reaches  its  destination.  Packets  of  the  same  con¬ 
versation  usually  follow  the  same  path.  However,  they  may  follow 
different  paths,  possibly  arriving  at  the  destination  out  of  order. 

In  this  case  they  would  have  to  be  reassembled  in  their  right  order 
before  being  delivered  to  the  ultimate  user.  Packet  switching  thus 
avoids  dedicating  the  network  facilities  to  users  during  silent  periods 
since  no  packets  are  generated  then.  In  hybrid  switching,  both  circuit- 
switched  and  packet-switched  traffics  exist  concurrently  on  the  same 
channel.  Each  channel  is  partitioned  into  two  subchannels;  a  circuit- 
switched  subchannel  and  a  packet-switched  subchannel.  Voice 
traffic  uses  the  circuit-switched  subchannel  while  data  traffic  uses 
the  packet-switched  subchannel.  To  increase  the  channel  utilization 
a  "movable  boundary"  feature  can  be  incorporated,  which  permits  the 
data  packets  to  use  any  residual  circuit-switched  capacity  that  may 
be  momentarily  available  due  to  voice  traffic  variations  [9]. 

A  well-designed  communications  network  ensures  that  the  traffic 
it  handles  does  not  exceed  the  capacity  of  the  network  devices  in  order 
that  congestion  may  not  build  up  and  reduce  the  throughput  of  the 
network.  In  the  extreme  case  congestion  may  lead  to  a  deadlock  in 
which  case  communication  becomes  impossible.  The  term  flow  control  is 
used  to  describe  any  mechanism  that  ensures  that  the  traffic  entering 
into  the  network  is  maintained  within  limits  compatible  with  the  capacity 
of  the  network  devices.  A  flow  control  protocol  should  perform  the 
following  functions  [16]: 


(a)  It  should  prevent  throughput  degradation  and  loss  of  efficiency 


due  to  overload. 

(b)  It  should  prevent  deadlock. 

(c)  It  should  ensure  a  fair  allocation  of  the  network  resources 
among  competing  users. 

(d)  It  should  ensure  that  the  rate  at  which  the  network  accepts 
data  is  matched  to  the  rate  at  which  it  can  transmit  them. 

Voice  and  data  traffics  make  different  and  conflicting  demands 
on  the  network:  They  show  different  tolerances  to  delay  and  errors. 

Voice  conversations  require  continuous  and  almost  real-time  delivery, 
they  are  very  sensitive  to  delay.  Data  conversations,  on  the  other 
hand,  are  generally  intolerant  of  errors  but  less  sensitive  to  delay. 

Data  conversations  need  to  be  very  reliably  delivered  to  their  ultimate 
users.  A  flow  control  protocol  that  is  to  be  applied  to  an  integrated 
voice  and  data  network  must,  therefore,  perform  two  additional 
important  functions,  namely: 

(a)  It  should  ensure  that  a  voice  conversation  is  transmitted 
with  essentially  constant  delay  while  trading  the  speech  quality  in 
response  to  nework  fluctuations. 

(b)  It  should  ensure  that  data  conversations  are  transmitted 
with  maximum  reliability  (that  is,  no  errors  or  lost  information)  with 
communications  delay  as  a  secondary  considerations. 

1 .3  Summary  of  Previous  Work 

The  concept  of  integrated  voice  and  data  communications  networks 
has  generated  much  interest.  It  has  been  motivated  by  the  need  to  exploit 
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recent  technological  advances  to  achieve  economies  of  scale:  sharing 
network  resources  more  efficiently  thereby  reducing  the  cost  per  user. 
Unfortunately  much  of  the  work  done  in  this  area  so  far  exists  as  mere 
heuristics  that  depend  on  simulation  for  evaluation. 

Coviello  and  Vena  [9]  propose  a  type  of  hybrid-switched  network 
called  the  slotted  envelope  network  (SENET)  which  provides  some  degree 
of  flexibility  for  a  varying  mix  of  traffic.  It  is  a  framed  TDM  scheme 
in  which  frames  are  divided  into  (V+D)  slots.  V  slots  are  allocated 
to  voice  conversations  on  a  circuit-switched  basis,  and  D  slots  are 
allocated  to  data  conversations  on  a  packet- switched  basis.  Voice  calls 
are  blocked  when  all  the  V  slots  reserved  for  voice  conversations  are 
busy.  A  voice  conversation  uses  its  slot  throughout  the  duration  of  the 
conversation.  In  order  to  increase  the  channel  utilization,  a  movable 
boundary  scheme  is  used  whereby  data  traffic  could  use  any  voice  slots 
that  may  be  momentarily  free  due  to  variations  in  the  number  of  active 
voice  conversations.  Forgie  and  Nemeth  [11]  consider  a  network  concept 
called  the  packetized  virtual  circuit  (PVC)  which  combines  selected 
features  from  the  packet  and  circuit  switching  technologies.  The  PVC 
technique  is  more  or  less  a  packet  switching  technique  in  which  messages 
are  packetized  and  travel  along  virtual  connections  to  their  destinations. 
The  authors  also  suggest  some  ways  to  deal  with  congestion  in  the  net¬ 
work,  including  discarding  some  voice  packets  when  queues  become  excessive. 
Bially  et  al.  [3]  propose  a  packet-switched  network  in  which  a  speech 
digitization  concept  called  "embedded  coding"  enables  the  network  to 
respond  dynamically  to  the  loading  conditions.  Here  speech  is  encoded 
into  different  packets  of  different  orders  of  importance.  If  all  the 


packets  are  delivered  to  the  sink,  the  result  would  be  a  high  quality 
speech.  When  the  network  loading  becomes  high  the  less  important  packets 
may  be  discarded.  The  delivered  packets  still  produce  usable  but  lower 
quality  speech.  Gitman  and  Frank  [17]  present  an  elaborate  economic 
analysis  of  the  integrated  voice  and  data  networks.  They  conclude,  on 
the  basis  of  a  number  of  assumptions,  that  even  when  no  optimization 
is  performed  in  selecting  the  voice  packet  length,  packet  switching  is 
still  superior  to  all  other  switching  techniques. 

1.4  Synopsis  of  the  Thesis 

The  hybrid-switched  network  is  basically  two  distinct  networks  in 
one,  if  the  movable  boundary  feature  is  not  incorporated.  Since  operational 
circuit-switched  and  packet-switched  networks  exist,  designing  a  pure 
hybrid-switched  integrated  voice  and  data  network  (i.e.  one  without 
movable  boundaries)  presents  no  new  challenges.  When  the  movable  boundary 
feature  is  incorporated,  the  difference  between  the  hybrid-switched  and 
the  packet-switched  network  becomes  a  matter  of  semantics  rather  than 
technology  as  the  operation  of  the  two  techniques  are  almost  identical, 
but  not  quite  anyway. 

We  design  a  packet-switched  integrated  network  that  achieves  a 
better  resource  allocation  than  the  pure  hybrid-switched  network.  We 
chart  a  new  course  by  formulating  the  integrated  network  problem  in  a 
way  that  lends  itself  to  an  analytic  and  algorithmic  solution.  Voice 
packets  have  a  non- preemptive  priority  over  data  packets.  The  reason 
for  practising  a  non- preemptive  priority  includes  the  fact  that  when  a 
data  packet  receiving  service  is  preempted  by  an  arriving  voice  packet. 
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the  part  of  the  data  packet  which  had  already  been  transmitted  before 
preemption  might  appear  as  errors  in  the  network.  We  could  avoid  this 
confusion  by  practising  a  non-preemptive  priority.  Also  it  does  not 
make  sense  to  use  network  resources  to  transmit  data  bits  that  would 
later  be  discarded. 

As  the  title  of  this  thesis  suggests,  we  are  primarily  interested 
in  the  flow  control  and  routing  issues  of  the  integrated  network.  We 
shall,  however,  discuss  some  protocols  for  congestion  avoidance  and 
control,  and  for  admitting  new  conversations  into  the  network. 

The  remainder  of  the  thesis  is  organized  as  follows.  In  Chapter  II 
we  present  the  formulation  of  the  network  model  and  the  problem  statement. 
In  Chapter  III  we  give  a  solution  to  the  problem  of  Chapter  II.  We  also 
formulate  a  joint  flow  control  and  routing  algorithm.  In  Chapter  IV  we 
consider  protocols  for  performing  miscellaneous  functions,  including 
congestion  avoidance  and  control,  and  the  admission  of  new  conversations 
into  the  network.  In  chapter  V  we  discuss  how  the  model  can  be  generalized 
for  a  multi-traffic  network  that  carries  traffic  with  different 
levels  of  delay  sensitivity.  In  Appendices  A  and  B  we  give  proofs  of 
theorems. 


CHAPTER  II 


MODEL  FORMULATION  AND  PROBLEM  STATEMENT 


2.1 .  Introduction 

In  [12]  Gallager  formulates  a  quasi-static  routing  strategy  for  a 
data  communication  network.  He  considers  a  network  with  slowly  varying 
average  input  rate  for  each  data  conversation,  and  develops  an  algorithm 
for  finding  minimal  delay  routing  using  distributed  computation.  Golestaani 
[19]  extends  Gallager's  work  by  considering  a  network  with  controllable 
input  rates.  These  rates  are  adjusted  periodically  on  the  basis  of  the 
information  on  the  network  utilization.  Also  these  rates  are  controlled 
in  a  manner  that  ensures  that  all  users  are  fairly  treated.  The  algorithm 
developed  in  [19]  implements  this  fairness  by  attempting  to  establish  an 
optimal  tradeoff  between  each  user's  satisfaction  and  the  network  loading 
condition.  A  user's  satisfaction  is  related  to  the  rate  at  which  he  is 
permitted  to  transmit:  the  closer  the  permitted  rate  is  to  his  desired 
rate,  the  more  satisfied  the  user  is. 

The  work  done  here  is  a  further  generalization  of  the  work  of 
Gallager  [12],  with  applications  to  the  integrated  voice  and  data  network. 

As  in  [19]  we  consider  an  integrated  network  with  controllable  voice 
and  data  input  rates.  The  voice  input  rate  is  controlled  by  adjusting 
the  voice  packet  lengths.  The  length  of  the  voice  packets  determines  the 
quality  of  the  voice  conversation:  the  longer  the  packets,  the  higher  is 
the  quality  of  the  conversation.  The  voice  packet  lengths  and  the  data 
input  rates  are  also  controlled  in  a  manner  that  ensures  that  all  users 


are  fairly  treated.  We  effect  this  fairness  by  developing  an  algorithm 
that  attempts  to  establish  an  optimal  tradeoff  between  the  network 
conditions  and  all  users’  satisfaction.  For  voice  conversations,  the 
longer  the  packets  (and  hence  the  higher  the  speech  quality)  the  more 
satisfied  the  user  is.  And  for  data  conversations,  the  closer  a  user's 
allowed  rate  is  to  his  desired  rate,  the  more  satisfied  he  is.  We  assume 
that  a  voice  digitizer  that  can  set  the  voice  packets  to  any  desired  lengths 
can  be  found.  We  do  not  concern  ourselves  with  what  happens  to  the 
speech  at  the  receiver  if  the  voice  packet  length  is  continuously  varied. 
This  is  an  aspect  of  signal  processing  and  has  not  yet  been  resolved; 
it  is  a  topic  for  future  research. 

The  work  done  here  is  not  the  first  to  consider  a  scheme  where 
voice  packet  lengths  are  varied  in  response  to  network  conditions.  The 
work  of  Bially  et  al .[3]  is  similar  to  ours.  The  way  flow  control  is 
effected  in  [3]  is  equivalent  to  varying  the  voice  packet  length  inside 
the  network,  and  their  scheme  is  more  dynamic  than  ours.  In  this  thesis 
we  apply  flow  control  at  the  input  to  the  network  by  making  changes  in 
voice  packet  lengths  before  sending  the  packets  into  the  network.  This 
would  then  reduce  the  need  to  drop  any  segments  of  these  packets  in  the 
network  as  is  equivalently  done  in  [3].  We  shall  discuss  the  similarity 
betweeen  the  two  models  more  in  Chapter  IV. 

2.2  The  Model 
2.2.1  Definitions 

Consider  a  store-and- forward  network  with  N  nodes  represented  by 
i  *  1,2,...,N.  A  link  that  goes  from  node  k  to  node  l  is  represented 
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by  the  symbol  (k,&),  We  assume  that  link  ( k ,Jt)  is  a  directed  link  and 
'is  thus  different  from  link  (t, k).  Since  most  voice  and  data  links  in 
use  are  full  duplex,  we  assume  that  if  link  (k,Jt)  exists  then  link  (t,k) 
also  exists.  Let 


f’kU) 


=  the  expected  voice  traffic  on  link  (i,k),  in 
bits/second,  destined  for  node  j 
=  the  expected  data  traffic  on  link  (i,k),  in 
bits/second,  destined  for  node  j 
=  the  expected  length  of  voice  packets,  in  bits, 
belonging  to  conversations  which  enter  the  net¬ 
work  at  node  i  and  are  destined  for  node  j 
*  the  expected  data  input  rate  at  node  i  ,  in 
bits/second,  of  conversations  destined  for  node  j 
=  the  rate  at  which  voice  packets  are  emitted  by 
the  voice  digitizer  during  a  talkspurt 


11 

5  fik^ 


—  the  aggregate  voice  traffic,  in  bits/second,  on 
link  (i,k) 

N 

=  Z  ft  (J) 

j=l 

=  the  aggregate  data  traffic,  in  bits/second, 
on  link  (i,k) 

s  F^k  +  ,  is  the  total  traffic  on  link  (i,k) 
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C.^  =  the  capacity  of  link  ( i , k)  in  bits/second 

0(i)  =  the  set  of  nodes  with  links  going  out  from  node  i 

I (i )  =  the  set  of  nodes  with  links  going  into  node  i 

m--  *  the  number  of  data  conversations  entering  the  network 

'  J 

at  node  i  and  destined  for  node  j 
n..  s  the  number  of  "off-hook"  speakers  with  voice  conversations 

■  J 

entering  the  network  at  node  i  and  destined  for  node  j 
n. .  is  made  up  of  two  parts: 

*  J 

nl.  =  the  number  of  n..  in  talkspurt 

1  J  1  J 

n^.  =  the  number  of  n..  in  silence, 

ij  ij 

When  a  speaker  is  off-hook  (i.e.  in  the  conversational  mode)  he 
alternates  randomly  between  the  talkspurt  and  the  silence  modes.  For 
a  large  value  of  n..,  the  number  of  speakers  in  talkspurt  can  be  modelled 

■  J 

by  a  Poisson  arrival  process  with  rate  9.  For  ease  of  analysis  we 
assume  that  the  holding  times  of  talkspurts  ar$  exponentially  distributed 
with  mean  <J>"\  Then  the  speaker  activity  for  the  (i,j)  off-hook  conver¬ 
sations  can  be  modelled  by  the  Markov  chain  shown  in  Fig.  2.1. 


Figure  2.1  Speaker  Activity  Model  for  (i,j)  Conversations 


Thus  the  activity  of  each  off-hook  speaker  can  be  modelled  by  the  two- 
state  Markov  chain  shown  in  Fig.  2.2- 
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T  =  Talkspurt  State 
S  =  Silence  State 


Figure  2.2  A  Markov  Chain  Model  of  A  Speaker 


The  steady-state  probability  that  the  speaker  is  in  talkspurt  is: 


*  '  fTT  <2-'> 

Therefore,  if  there  are  n. .  off-hook  speakers,  the  steady-state 

*  J 

probability  that  k  of  them  are  in  talkspurt  is: 


Pk 


k 

Y 


(1 


(2.2) 


One  problem  in  the  modelling  of  a  voice  network  is  that  for  any  given 
number  of  off-hook  speakers  the  number  of  speakers  in  talkspurt  changes 
too  fast  to  be  tracked  by  any  reasonable  algorithm  that  makes  use  of 
global  information.  We  are  then  left  with  the  choice  of  estimating  nl., 

*  J 

given  n^j.  In  this  work  we  initially  estimate  nlj  by  the  mean  value  of  k; 


-15- 


2.2.2  Voice  Packet  Delay  Variations 

Traditional  circuit  switching  systems  exhibit  fixed  delay  to 
voice  traffic.  Packet  switching  systems,  on  the  other  hand,  exhibit 
variable  delay  due  to  the  possible  queueing  of  the  voice  packets  at  the 
different  nodes  on  the  route  from  the  source  to  the  destination.  It  is 
assumed  that  when  a  speaker  is  in  talkspurt  the  voice  digitizer  emits 
voice  packets  at  a  constant  rate  of  8  packets  per  second  [3].  These 
packets  experience  different  queueing  delays  along  their  route.  If  each 
packet  is  delivered  to  the  ultimate  user  (or  sink)  as  soon  as  it  arrives 
at  the  destination  node,  the  received  message  would  contain  many  uneven 
and  annoying  gaps.  Therefore,  it  is  necessary  that  we  install  a  smoothing 
buffer  at  the  destination  node,  temporarily  store  the  arriving  voice 
packets,  and  finally  release  them  to  the  sink  at  the  same  rate,  8 
packets/second,  at  which  the  voice  digitizer  sent  them  into  the  network. 
In  Chapter  IV  we  shall  consider  the  case  where  the  rate  at  which  the 
packets  are  released  to  the  sink  is  different  from  8. 

A  schematic  diagram  of  what  happens  to  the  voice  packets  as  they 
travel  along  the  route  from  source  node  i  to  destination  node  j  is 
shown  in  Fig.  2.3.  Figure  2.4  is  a  diagram  of  the  possible  path  delays 
of  five  packets  of  a  conversation.  The  first  packet  is  transmitted  at 
time  zero;  the  second  packet  at  time  1/8;  and  in  general,  the  kt*1  packet 
is  transmitted  at  time  k  >_  1.  There  are  three  links  on  the  path. 


From  other  nodes 


From  the  diagram  we  observe  that  the  first  packet  suffers  a  delay  dQ, 
since  it  arrives  at  time  d  .  This  packet  is  not  delivered  immediately 
to  the  sink  but  is  held  back  in  the  buffer  for  some  time  t  -  drt,  and 
delivered  at  time  t  .  The  reason  for  this  action  will  be  clear  shortly. 
The  second  packet  arrives  at  time  d^ ;  it  is  held  back  in  the  buffer  and 
released  at  time  tq  +  1/8.  The  third  packet  arrives  at  time 
d^  =  rQ  +  2/3,  just  in  time  to  be  delivered  to  the  sink.  The  fourth 
packet  arrives  later  than  it  is  required  to  be  delivered  to  the  sink 
and  so  is  discarded.  A  protocol  could  be  formulated  for  generating 
fictitious  packets  which  are  to  be  released  to  the  sink  when  the  actual 
packets  arrive  late.  However,  this  is  a  signal  processing  function  which 

fU 

we  will  not  be  concerned  with  here.  In  general,  if  the  k  packet 

k-1 

arrives  earlier  than  t  +  ,  k  >  1,  it  is  stored  in  the  buffer 

o  S  — 

k-1 

and  released  to  the  sink  at  time  t  +  .  If  it  arrives  later  than 

o  3 

this  time,  it  is  thrown  away. 

The  reason  for  imparting  the  extra  delay  on  the  first  packet  is 
now  obvious:  to  help  reduce  the  likelihood  of  discarding  other  packets 
due  to  their  late  arrival.  Obviously  choosing  the  extra  delay  tq  -  dQ 
to  be  very  large  ensures  that,  with  a  probability  close  to  one,  all 
packets  arrive  in  time  to  be  delivered  to  the  sink.  However,  this  is 
not  the  solution  we  seek.  We  have  noted  that  voice  conversations  require 
almost  real-time  delivery;  therefore,  the  recommended  choice  of  the  extra 
delay  is  tq  -  dQ  £  1/8. 

One  of  our  objectives  in  the  above  scheme  would  be  to  reduce  the 
path  delay  of  each  voice  packet.  This  would  in  turn  reduce  the  likelihood 
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of  the  packet  being  discarded.  Another  objective  would  be  to  increase 
the  qua! it-  -f  each  voice  conversation.  A  realization  of  these  objectives 
resides  in  making  tradeoffs  between  speech  quality  and  delay.  Short 
packets  would  travel  faster  than  long  packets  for  two  reasons.  First, 
there  is  less  transmission  delay  per  packet  on  each  link.  Secondly, 
short  packets  generate  smaller  traffic  than  long  packets,  and  so  there 
is  reduction  in  congestion.  Thus  the  use  of  short  packets  would  solve 
the  delay  problem-.  However,  short  packets  imply  poorer  quality  speech. 
Therefore,  we  need  to  use  long  packets  in  order  to  generate  higher 
quality  speech;  but  there  is  the  attendant  possibility  of  incurring 
larger  packet  delay.  Any  attempt  to  gain  in  quality,  therefore,  results 
in  poorer  delay  performance;  and  any  attempt  to  improve  the  delay  performance 
degrades  the  quality.  We  need  then  to  specify  what  would  be  the  point  of 
optimal  tradeoff  between  quality  and  delay  such  that  once  we  attain  it, 
we  would  not  be  tempted  to  improve  the  quality  or  the  delay  performance 
any  further.  This  suggests  that  we  formulate  an  optimization  problem 
the  solution  of  which  would  generate  the  optimal  packet  lengths  (and 
hence  quality)  and  the  optimal  delay.  Our  task  then  boils  down  to 
choosing  the  appropriate  objective  function  for  the  problem.  Since  we 
are  dealing  with  an  integrated  voice  and  data  network,  we  must  ensure 
that  the  solution  to  the  problem  we  shall  formulate  also  provides  the 
optimal  data  input  rates  and  the  optimal  data  packet  delay.  In  the 
subsequent  subsections  we  discuss  the  different  aspects  of  the  problem, 
and  find  appropriate  cost  functions  that  will  constitute  the  objective 
function. 
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2.2.3  Cost  of  Link  Congestion 

Let  dj^  =  the  delay  experienced  by  the  packet  of  a  talkspurt  on 

link  (k,$)  including  procesing  and  queueing  delays  at  node  k 

R  =  the  route  of  the  voice  conversation 

Then  dn  *  ^  dj^  ,  is  the  path  delay  of  the  nth  packet  of  the 

( k  ,Jt)  eR 

talkspurt. 

In  the  discussion  above  the  delay  experienced  by  the  first  packet 
of  a  talkspurt,  d^ ,  is  crucial  in  determining  whether  or  not  a  subsequent 
packet  is  discarded.  If  the  first  packet  experiences  little  delay  while 
other  packets  experience  larger  delays,  then  many  of  these  packets 
would  be  discarded.  On  the  other  hand,  if  the  first  packet  experiences 
a  large  delay  while  other  packets  experience  less  delay,  then  not  many 
of  these  packets  would  be  discarded.  Thus  the  probability  that  the  n 
packet  of  a  talkspurt  would  be  discarded  depends  not  only  on  dn  but  also 
on  d\  This  dependence  on  d^  is  not  a  very  desirable  feature.  In  order 
to  reduce  this  dependence  (it  cannot  be  totally  eliminated!),  we  will 
attempt  to  make  all  packets  of  the  same  talkspurt  experience  similar 
delays  on  each  link  on  their  path.  This  would  then  translate  to 
their  experiencing  similar  path  delays,  and  hence  reducing  the  chances 
of  their  being  discarded. 

We  hope  to  realize  the  above  objective  as  follows.  We  consider  a 
cost  to  be  associated  with  the  link  delay  of  each  voice  packet.  Similarly, 
we  consider  a  cost  to  be  associated  with  the  link  delay  of  each  data  packet. 
Since  the  delay  requirements  of  voice  packets  are  more  stringent  than 


those  of  data  packets,  we  grant  the  voice  packets  a  non-preempti ve 
priority  over  the  data  packets,  Packets  wait  on  a  first-come  first- 
served  basis  within  their  priorities.  We  assume  that  the  queue  at  each 
output  link  can  be  modelled  by  an  M/G/l  queue  with  priorities  [22]. 

Consider  a  tagged  packet  that  arrives  at  a  node  for  service. 

Let 

Ay  =  the  Poisson  rate  at  which  voice  packets  arrive  at  the 
output  queue  of  the  node,  in  packets/second 

=  the  Poisson  rate  at  which  data  packets  arrive  at  the 
queue,  in  packets/second 

ny  =  the  number  of  voice  packets  already  in  queue  when  the 
tagged  packet  arrived 

n^  =  the  number  of  data  packets  already  in  queue  when  the 
tagged  packet  arrived 

I 

ny  =  the  number  of  voice  packets  which  arrive  while  the  tagged 
packet  is  in  queue. 

Note  that  if  the  tagged  packet  is  a  voice  packet  these  ny  packets  have 
no  effect  on  its  waiting  time. 

Tr  =  the  remaining  time  to  complete  the  service  of  the  packet 
in  service  when  the  tagged  packet  arrived 
=  the  time  to  service  all  the  ny  yoice  packets 

•  i 

Tl  =  the  time  to  service  all  the  ny  voice  packets 

Tg  =  the  time  to  service  all  the  n^  data  packets. 


Then  if  the  tagged  packet  is  a  voice  packet,  its  waiting  time  is: 


Taking  expectations  we  obtain 

Wy  =  +  W  -|  (2.4) 

where  W.  =  E(T^).  If  the  time  to  service  one  voice  packet  is  ry,  then 

Wt  ■  ECny)  E(tv)  =  E(ny)  ry  . 


By  Little's  formula  [24],  E(nd)  =  Xd  W^.  We  cannot  apply  Little's  formula 
directly  to  E(ny),  but  we  can  establish  its  value  by  the  following  argument. 
Since  we  know  that  a  data  packet  waits  for  an  average  time  Wd  and  that  voice 
packets  arrive  at  an  average  rate  Xv  ,  clearly  the  average  number  of 
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voice  packets  which  arrive  over  the  time  interval  Wd  is  E(ny)  =  Xy  Wd- 
A  more  rigor  ,  derivation  of  this  formula  can  be  found  in  [22], 
Therefore,  from  [2. 6) 


S  Wj.  +  Xy  Wy  Ty  +  Ay  Wd  Ty  + 


or  _  *r  +  Py  My  _  *r+  Pv  Wv 

Wd  ~  T  -  Py  -  Pd  “  1  -  p 

A  _ 

where  pd  =  Xd  Td  is  the  link  data  utilization  factor 
A 

p  =  Py  +  pd  <  1,  is  the  link  aggregate  utilization  factor. 

Substituting  for  Wy  we  obtain 


We  may  rewrite  (2.5)  and  (2.7)  as  follows: 

wrc  wrc 

u  s  — _  =  — - 

v  C-o  C  C-F 

^v  v 


(2.8) 


where  Fy  =  the  link  voice  traffic, in  bits/second, 

F  *  the  aggregate  link  traffic,  in  bits/second, 

C  a  the  link  capacity,  in  bits/second. 

For  completeness,  the  expected  remaining  time  to  completion  of  service 
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of  the  packet  in  service  when  the  tagged  packet  arrived  is  given  by 

[22]: 

Wr  -  5"  CXV  E(t^)  +  Xd  (2.10) 

The  delay  of  a  voice  packet  at  a  link  (&,k)  is  made  up  of  the 
waiting  time  and  the  service  time  (including  queueing  and  processing 
delays  at  node  £).  The  latter  is  usually  independent  of  the  network 
status.  Hence  at  each  node  we  bother  only  about  the  waiting  time  at  the 
output  link.  From  (2.8)  we  observe  that  the  expected  waiting  time  of 
a  voice  packet  is  a  function  of  the  link  voice  traffic  F  .  This 
expected  waiting  time  is  also  affected  by  the  amount  of  time  the  voice 
packet  has  to  wait  for  a  data  packet  already  in  service  when  the  voice 
packet  arrived;  but  this  effect  is  rather  small.  Because  of  these 
facts  we  define  a  cost  function  B^fY^)  as  follows:  B^F^)  is  the 
cost  of  limiting  the  voice  traffic  on  link  (i,k)  to  fY^.  For  mathematical 
tractability  we  assume  that  B^fY^)  is  a  convex  increasing  and  twice 
differentiable  function  of  fY^,  with  a  typical  plot  as  shown  in  Fig. 2. 5. 


The  reason  for  the  cutoff  at  is  the  following. 

From  (2.8)  we  observe  that  as  Fy  approaches  C,  the  waiting  time  becomes 
large.  Since  voice  packet  delay  is  a  critical  issue,  we  limit  the  voice 
traffic  on  each  link  to  some  predetermined  fraction,  ,  of  the  link 
capacity.  This  flow  restriction,  coupled  with  the  fact  that  voice 
traffic  has  a  non-preemptive  priority  over  data  traffic,  would  enable 
the  voice  packets  to  experience  less  delay  than  would  be  the  case  if 
the  restriction  were  removed.  That  is,  if  we  hold  the  voice  traffic 
below  a-|  C^,  the  link  delay  would  be  small. 

We  could  define  a  cost  function  for  the  data  traffic  in  a  similar 
manner.  However,  we  can  exploit  some  characteristics  of  voice  traffic  in 
defining  the  cost  function  for  the  data  traffic.  First,  observe  from  (2.9) 
that  the  expected  waiting  time  of  a  data  packet  at  the  output  link  is  a 
function  of  both  Fy  and  F.  One  major  attraction  of  the  integrated  voice 
and  data  network  is  the  anticipated  possibility  of  exploiting  the  on-off 
characteristics  of  voice  traffic  to  transmit  more  data  when  voice  traffic 
is  low.  A  reasonable  cost  function  for  data  traffic  should,  therefore, 
recognize  this  fact:  when  fY^  =  0  the  data  packets  should  be  able  to 
use  link  (i,k)  as  if  the  network  were  an  all-data  network.  Since  data 
packet  delay  requirements  are  not  as  stringent  as  those  for  voice  packets, 
we  would  permit  a  higher  cutoff  point  for  each  link  flow.  A  cutoff  is 
necessary  to  account  for  the  finite  buffer  spaces  available  at  the 
ojtput  links. 

That  data  packet  delay  is  more  tolerable  than  voice  packet  delay 
means  that  the  marginal  cost  to  a  link  of  a  voice  packet  should  be  at 
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least  as  great  as  that  of  a  data  packet.  That  is,  the  cost  of  congestion 
for  an  additional  voice  packet  is  greater  than  that  for  an  additional 
data  packet.  When  a  given  voice  packet  arrives  at  a  node,  it  holds  up 
the  service  of  all  data  packets  in  queue  on  its  arrival.  Its  delay  is 
only  affected  by  the  voice  packets  which  are  in  queue  on  its  arrival. 

Data  packets  affect  the  voice  packet's  delay  through  the  possibility 
that  a  data  packet  was  receiving  service  when  the  voice  packet  arrived. 

As  we  said  earlier,  this  effect  is  rather  small.  When  a  given  data 
packet  arrives  at  a  node,  its  delay  is  affected  not  only  by  all  the 
packets  (both  voice  and  data)  in  queue  on  its  arrival  but  also  by  all 
voice  packets  which  arrive  while  it  is  in  queue.  We  then  need  a  model 
that  recognizes  this  difference  in  the  marginal  costs  for  voice  and 
data  conversations. 

Now  consider  the  cost  function  k(Fi k) defined  as  follows:  Gik(Fik)is 
the  cost  of  limiting  the  aggregate  traffic  on  link  (i,k)  to  F^.  For 
reason  of  mathematical  tractability  we  assume  that  G-k(F^ k)  is  a  convex 
increasing  and  twice  differentiable  function  of  F^,  with  the  typical 


Let  the  composite  function  F^)  be  defined  as  follows; 


°1k<F1k*  F?k>  “  Bik(F1k)  +  G1k(Fik) 


(.2.11) 


That  the  composite  function  D-k  captures  the  essence  of  our  design 
objectives  for  congestion  as  listed  above  can  be  seen  from  the  following 
analysis.  When  *  0,  then  from  the  way  we  have  defined  B^F^)  we 
have  that  Bik(0)  *  0  and 

Gik<-FikiFTk  =  °>  •  sik(F?k> 

Thus 

°ik(Fik-°!  *  6ik<F?k>  <2-12> 

This  is  precisely  the  cost  function  defined  by  Golestaani  [19]  for  an 
all-data  network.  Also  the  derivatives 


3D 


ik 


»F« 


*  Bik  <Flk>  +  Sk  ‘Fik> 


a°,k 

’  Glk  <Fik> 


[where  Bik  (F^)  *  dBik  ,  etc.] 

dFIk 

are  the  marginal  costs  to  link  (i,k)  of  voice  and  data  conversations 
respectively.  As  can  be  seen  above, the  marginal  cost  of  voice  to  link 
(i,k)  is  greater  than  that  of  data,  as  we  hoped. 


In  suninary  then,  our  argument  is  that  we  need  two  cost  functions 
that  relate  to  delay  (or  congestion)  in  the  integrated  network.  One  cost 
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function  deals  with  the  issue  of  ensuring  that  we  Toad  the  network  with 
voice  traffic  to  such  a  level  that  would  not  make  the  voice  packet  delay 
excessive.  The  other  cost  function  deals  with  the  mix  of  voice  and  data 
traffics.  This  latter  function  allows  us  to  exploit  the  on-off  character¬ 
istics  of  voice  traffic  to  transmit  more  data  when  voice  traffic  is  low. 

2.2.4  Voice  Quality  Limitation  Cost 

Having  considered  the  issue  of  congestion  in  the  network  we  now 
deal  with  the  issue  of  fairness  to  all  voice  conversations.  We  wish 
to  design  a  network  that  operates  as  follows.  When  a  voice  conversation 
is  established  the  network  is  bound  to  accept  all  packets  of  that  conversa 
tion  as  long  as  the  speaker  is  off-hook.  However,  the  quality  of  that 
conversation  is  not  guaranteed  to  be  good  at  the  sink.  The  quality  of 
different  talkspurts  of  the  same  conversations  may  even  be  different, 
depending  on  the  status  of  the  network  when  a  talkspurt  is  generated. 

As  we  have  mentioned  earlier,  the  quality  of  any  voice  conversation  is 
related  to  the  length  of  the  packets  of  that  conversation:  the  longer  the 
packets,  the  higher  the  quality  of  the  conversation. 

In  order  to  effect  fairness  to  all  voice  conversations,  then,  we 
make  it  increasingly  costly  to  degrade  the  quality  of  any  conversation 
any  further  through  assigning  shorter  and  shorter  packets  to  that 
conversation.  Let  a  voice  conversation  that  enters  the  network  at 
node  1  and  is  destined  for  node  j  be  denoted  by  voice  conversation 
(i,j).  For  each  voice  conversation  (i,j)  we  define  a  nominal  packet 
length  a..;  that  is,  the  maximum  length  we  can  assign  to  packets  of 

*  J 

conversation  (i,j)  is  a...  The  reasons  for  this  packet  length  limitation 

•  yJ 
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are  as  follows.  First,  above  a  certain  packet  length  there  is  no 
appreciable  improvement  in  the  yoice  quality.  Therefore,  it  would  not 
make  sense  to  assign  lengths  that  would  generate  more  traffic  without 
any  extra  gain  in  quality.  Secondly,  the  voice  digitizer  is  limited 
in  the  number  of  bits  it  can  generate  over  a  given  time  interval.  For 
each  voice  conversation  (i,j)  we  define  a  cost  function  V..(i..)  as 

'  J  *  J 

follows:  V..(i..)  is  the  cost  of  restricting  the  packet  length  of  voice 

I  J  I  J 

conversation  (i,j)  to'  £...  For  mathematical  tractability  we  assume 

•  v 

that  is  a  convex  non-increasing  and  twice  differentiable  function 

of  with  the  typical  plot  as  shown  in  Fig.  2.7.  For  notational 
*  J 

convenience  we  assume  that  all  voice  conversations  (i,j)  have  the  same 
cost  function  V^U^). 


Figure  2.7  Quality  Limitation  Cost  of  A  Voice  Conversation  (i,j) 
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2.2.5  Data  Rate  Limitation  Cost 

For  data  conversations  we  let  the  allowable  input  rates  depend 
on  the  network  loading  conditions:  when  the  network  is  heavily  loaded 
we  make  the  input  rates  small,  and  when  the  network  is  lightly  loaded 
we  permit  higher  input  rates.  We  thus  assume  that  the  lengths  of  data 
packets  are  independent  of  the  network  control.  In  order  to  effect 
fairness  to  all  data  conversations  we  make  it  increasingly  costly  to 
cut  back  any  conversation  further  through  assigning  it  a  smaller  and 
smaller  input  rate.  In  this  regard  we  define  a  cost  function  E..(r..) 
as  follows:  E..(r.  .)  is  the  cost  of  restricting  the  input  rate  of 

'  J  *  \J 

data  conversation  (i,j)  to  r  For  reason  of  mathematical  tractability 

*  J 

we  assume  that  it  is  a  convex  non- increasing  and  twice  differentiable 
function  of  r_  ,  with  the  typical  plot  as  shown  in  Fig.  2.8.  For 
simplicity  of  notation,  we  take  E..(r .  . )  to  be  the  same  for  all  data 

'  J  '  J 

conversations  (i,j).  r..  is  some  desired  rate  at  which  data  conversation 

'  J 

(i,j)  would  be  transmitted  if  we  did  not  exercise  any  control  over  it. 


2.3  The  Problem  Statement 

Our  model  formulation  is  now  complete.  The  objective  function 


of  our  problem  is 

a  sum  of  four  terms: 

CD 

bt(fv) 

- 

5  Bik(F"k) 

cm 

Gt(.F) 

- 

£  Gik<Fiki 

(in) 

vt(D 

- 

(iv) 

ET(r) 

= 

^  Dj  Eij(rij) 

where  Yn..  is,  as  we  stated  earlier,  our  estimate  of  nl.,  the  number 
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of  off-hook  speakers  in  talkspurt;  and  m..  is  the  number  of  data 

i  <J 

conversations  (i,j).  We  can  then  formulate  the  problem  as  follows: 


Minimize  J 
subject  to 


Bt(Fv)  +  Gt(F)  +  Vt(M  +  Ey(r) 


teO(i) 
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(2.13] 

(2.14] 

(2.15] 
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Constraints  (2.14)  and  (2.15)  are  the  so-called  continuity  (or  flow 
conservation)  equations.  They  state  that  the  total  traffic  coming  into 
node  i  and  destined  for  node  j  is  equal  to  the  total  traffic  going 
out  of  node  i  and  destined  for  node  j.  Constraint  (2.16)  states  that 
flows  are  non-negative.  In  the  next  chapter  we  present  a  solution  to 
this  problem  and  develop  a  joint  flow  control  and  routing  algorithm  that 
sets  the  voice  packet  lengths  £..  and  data  input  rates  r.  •  to  values 

I  J  I  J 

compatible  with  the  network  status,  and  for  routing  the  voice  and  data 
traffics. 


CHAPTER  III 


SOLUTION  TO  THE  PROBLEM 

3.1  Introduction 

In  Chapter  II  we  formulated  a  model  of  the  integrated  voice  and 
data  network  as  a  convex  optimization  problem.  In  this  chapter  we 
present  a  solution  to  this  problem.  We  construct  a  joint  flow  control 
and  routing  algorithm  that  sets  the  voice  packet  lengths  and  the  data 
input  rates,  and  routes  the  voice  and  data  traffics,  in  a  manner  that 
attempts  to  satisfy  the  optimality  conditions  we  shall  derive. 


3.2  The  solution 

Theorem  3.1  Let  u  =  ,  r _ ,  f Y k ( j ) ,  f?k(j)}  be  a  feasible  point 

of  (2.14)  through  (2.16).  Then  u  minimizes  (2.13)  if  and  only  if  there 
exist  two  sets  of  numbers  X  =  (X. .}  and  y  =  {y..}  »  with  X..  =  0  and 
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y..  =  0,  such  that  the  following  Kuhn-Tucker  conditions  are  satisfied. 
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(3.4) 

where  denotes 

the  derivative  of 

Bik^Fik^’  etc-  The  Proof  of  this 

theorem  is  given  in  Appendix  A. 


3.2.1  Definitions.  A  voice  route  rY.  between  nodes  i  and  j  is  a  set 

’  J 

of  links  ((i,k).  (k,i) . . » (m, j ) }  which  connects  nodes  i  and  j  and 
along  which  voice  traffic  flows  from  node  i  to  node  j.  We  define  a 
data  route  R?j  in  a  similar  manner. 

We  can  interpret  the  quantity  +  Gi k(Fi k)  as  the  marginal 

cost  of  voice  traffic  on  link  (i,k),  or  the  "voice  length"  of  link  (i,k). 
Similarly,  the  quantity  G 1 k ( F^ k)  can  be  interpreted  as  the  marginal  cost 

of  data  traffic  on  link  (i.k)  or  the  "data  length"  of  link  (i,k).  Since 
Xjj  =  0  and  pjj  =  0,  we  can  solve  for  X^j  and  u..  recursively  to  obtain 
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Then  we  can  interpret  X--  as  the  marginal  voice  cost  of  congestion  on 
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a  path  of  optimal  flow,  or  the  voice  length  of  path  R^.  when  flow  is 
optimal.  Similarly,  we  can  interpret  p. .  as  the  marginal  data  cost  of 
congestion  on  a  path  R?.  of  optimal  flow,  or  the  data  length  of  path 

v 

R^.  when  flow  is  optimal.  Equations  (3.1)  and  (3.2)  then  state  that  all 
traffic  flows  along  paths  of  minimum  marginal  cost,  i.e.  "shortest"  paths. 
In  (3.5)  and  (3,6)  and  p . .  are  defined  as  properties  of  the  optimal 

i  J  l  J 

flow  and  are  difficult  to  find  without  solving  the  optimization  problem. 

We  shall  need  a  form  of  X. ■  and  p. .  in  our  joint  flow  control  and  routing 
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algorithm.  Therefore,  we  shall  redefine  X- ■  and  p..  in  terms  of 
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measurable  quantities,  as  functions  of  an  arbitrary  flow  as  follows: 


X  . . 
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""  E  #  ^  <Fu>  +  sn<Fn>3 
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(3.7) 


(3.8) 


When  flow  is  optimal,  the  X..  obtained  from  (3.5)  is  equal  to  that 
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obtained  from  (3.7);  and  theP..'s  obtained  from  (3.6)  and  (3.8)  are  equal. 
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3.2.2  The  Priority  Functions  [19] 

We  define  the  voice  priority  function,  p^U^),  for  voice 
conversation  (i,j)  as  follows: 


-  r  <l«) 


Similarly,  we  define  the  data  priority  function,  q^tr.^),  for  data 


(3.9) 


conversation  (i ,  j )  as 


(3,10) 


Then  we  can  restate  the  optimality  conditions  (3.3)  and  (.3.4)  as  follows: 
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Recall  that  in  section  2.2.2  we  stated  that  the  solution  to  the  integrated 
network  problem  resides  in  finding  an  optimum  compromise  between  speech 
quality  and  voice  delay,  and  between  data  input  rate  and  network  congestion. 
The  priority  functions  have  an  interesting  interpretation  in  this  regard. 
First,  the  priority  functions  p. .,(£..)  are  the  marginal  gain  in  voice 
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quality  for  an  additional  voice  packet  length  allocation.  What  (3.11) 
states  is  that  optimality  occurs  when  the  marginal  gain  in  voice  quality 
is  as  close  as  possible  to  the  marginal  voice  cost  of  congestion  A. ., 

*  v 

★  .  # 
subject  to  0<  t.,  <  a...  The  equilibrium  point  (or  point  of  optimal 
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tradeoff ). for  any  voice  conversation  is  attained  when  its  marginal  gain 
in  quality  equals  its  marginal  cost  of  congestion.  At  any  other  point 
the  marginal  cost  of  congestion  of  that  conversation  will  be  either  greater 
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than  or  less  than  its  marginal  gain  in  quality.  This  means  that  any 
improvement  in  quality  at  the  non-optimal  point  would  be  obtained  at 
the  expense  of  voice  packet  delay,  and  vice  versa.  The  same  explanation 
holds  for  data  conversations.  A  desirable  flow  control  algorithm  would 
then  be  one  that  attempts  to  equalize  these  two  sets  of  marginal  values 
for  each  conversation.  In  section  3.2.3  we  shall  exploit  the  above 
fact  in  setting  the  voice  packet  lengths  and  the  data  input  rates.  In 
section  3.4  we  shall  consider  a  class  of  priority  functions. 

3.2.3  A  Joint  Flow  Control  and  Routing  Algorithm 

The  joint  flow  control  and  routing  algorithm  we  propose  performs 
four  different  network  functions:  it  effects  voice  flow  control  by 
setting  the  voice  packet  lengths  to  values  appropriate  for  the  network 
conditions;  it  effects  data  flow  control  by  setting  data  input  rates  to 
appropriate  values;  and  it  routes  voice  and  data  traffics  along  the  shortest 
paths.  There  are  two  approaches  to  defining  the  algorithm.  We  call  these 
approaches  the  "all -at-once  protocol"  and  the  "cyclic  coordinate  protocol". 

In  the  all-at-once  protocol,  the  algorithm  performs  all  the  four 
network  functions  simultaneously  during  each  iteration.  In  the  cyclic 
coordinate  protocol,  the  algorithm  performs  the  four  network  functions  in 
a  step-by-step  manner,  changing  one  parameter  (or  coordinate)  at  a  time. 

Specifically,  after  an  update  of  the  network  condition  is  made,  voice 

/ 

flow  control  is  effected;  another  update  is  made  after  a  steady  state  has 
been  reached,  and  voice  routing  changes  are  made,  Then  after  another 
steady  state  is  reached,  another  update  is  made  and  data  flow  control 
is  effected.  Then  again  another  update  is  made  after  another  steady  state 
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is  reached,  and  data  routing  changes  are  made.  The  process  is  repeated 
starting  with  another  update  followed  by  voice  flow  control.  Thus  there 
are  four  updates  to  complete  one  iteration  of  the  algorithm.  The 
convergence  of  the  cyclic  coordinate  algorithm  appears  to  be  slower  than 
that  of  the  all-at-once  algorithm.  We  shall  consider  only  the  all-at- 
once  algorithm  here. 

Our  joint  flow  control  and  routing  algorithm  belongs  to  the  class 
of  algorithms  proposed  by  Bertsekas  [2J.  This  class  of  algorithms 
operates  in  the  space  of  path  flows  rather  than  link  flows  in  which  our 
problem  is  defined.  However,  the  correspondence  between  a  problem  defined 
in  the  space  of  link  flows  and  one  defined  in  the  space  of  path  flows 
can  be  easily  established  [2]. 

An  iteration  of  the  algorithm  starts  with  each  node  i  computing 
the  shortest  voice  and  data  distances  from  itself  to  every  node  j  for 
which  node  i's  conversations  are  destined.  We  shall  explain  how  distance 
computation  is  carried  out  later.  The  metric  for  distance  computation 
is  the  marginal  cost  of  congestion  on  each  link  sunmed  over  all  links 
on  the  path  from  node  i  to  node  j.  Generally  there  may  be  two  or  more 
paths  that  qualify  for  the  shortest  path  for  each  (i,j)  pair.  Voice 
packet  lengths  and  data  input  rates  are  increased  or  decreased  according 
to  whether  the  priority  functions  are  greater  than  or  less  than  the 
shortest  distances.  Routing  changes  are  made  as  follows.  The  traffic 
on  each  non-shortest  path  Is  decreased  by  an  amount  proportional  to  the 
difference  between  the  length  of  that  path  and  the  length  of  the  shortest 
path.  The  traffic  on  the  shortest  path  is  then  the  total  traffic  generated 
by  the  (i,j)  conversations  less  the  traffic  routed  along  the  other  paths. 


If  two  or  more  paths  qualify  for  the  shortest  path,  we  arbitrarily  choose 
one  of  them  to  perform  the  function  described  above. 


Let  the  state  of  node  i  at  iteration  n  be  denoted  by 
un  *  Ujj  ,  r1?.,  ffkU)) 


(3.13! 


ij  ’  i J 

Since  there  may  be  more  than  one  shortest  path  for  voice  and/or  data 
conversations  (i,j),  we  define  the  joint  flow  control  and  routing 
algorithm  H  as  follows.  Let 

U  3  the  set  of  feasible  points  of  (2.14)  through  (2.16). 

Then  the  algorithm  H  is  a  point-to-set  mapping:  it  maps  a  point  uel) 
into  a  set  of  points  contained  in  U.  Thus  for  un  dJ,  the  algorithm  yields 
H ( u°) ^  U  from  which  an  arbitrary  element  un+^  is  selected.  Thus  given  an 
initial  point  u°,  the  algorithm  generates  sequences  through  the  iteration 
[25]: 


u"+1  c  H(un) 


(3-14) 


Let  ,*» 


6  y  n. .  i". 
T  U 


the  total  expected  rate  of  all  the  voice  conversation 
entering  the  network  at  node  i  and  destined  for  node  j 
at  iteration  n. 


.dn 


’lo¬ 


rn 

ij  ij 


the  total  expected  rate  of  all  data  conversations 
entering  the  network  at  node  i  and  destined  for 
node  j  at  iceration  n 


For  each  a  =  (i,j)  eA  ,  denote  by  svn(a) ;  and 
for  each  a  =  (i,j)  eA^,  denote  sdj  by  sdn(a) 

Pa  =  the  set  of  directed  paths,  with  no  repeated  nodes,  originating 

a 

at  node  i  and  terminating  at  node  j 
Spn  *  the  total  voice  flow  for  a  e  Ay  along  path  pcPa  at  iteration  n 

Sdn  *  the  total  data  flow  for  a  eA^  along  path  pePa  at  iteration  n 

Then  the  relationship  between  the  link  flows  and  Fd£  and  the  path 
flows  is  given  by 


t-vn 

Fik s 

E 

aeAy 

E  V 

P ePa 

i,k)  s™ 

(3.15) 

Fdn 

Mk 

E 

aeAd 

E  y 

pePa 

1,k)  sf 

(3.16) 

where  the  incidence  term<5p(i,k)  is  defined  as  follows: 

(  1  if  link  (i ,k)ep 

5p(i,k)  =  | 

(o  otherwise 

Let  "5  -  Bik  (f?S) +  4  <F?k> 

=  the  voice  length  of  link  (i,k)  at  iteration  n 
dik  “  Gik  (Fik^ 

*  the  data  length  of  link  (i,k)  at  iteration  n. 

Then  for  each  a  e  Ay  and  pePa  we  define  the  voice  length  of  path  p  at 
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iteration  n  as 


jvn 


ivn 


03.17) 


P  (t,k)ep 

Similarly,  for  each  a  eA^  and  p  ePa  we  define  the  data  length  of  path  p 


at  iteration  n  as 
,dn 


U,k)ep  £k 


(3.18) 


Finally,  we  define  the  shortest  voice  and  data  distances  from  node  i  to 
node  j  at  iteration  n  as 


jvn, 


.dn, 


-  Min  d™ 

(3.19) 

pcPa  » 

=  Min  d*1 

(3.20) 

pcPa  " 

respectively.  Many  efficient  methods  exist  for  finding  shortest  paths  in 
a  network;  some  of  these  methods  are  given  in  [23].  Then  the  joint  flow 
control  and  routing  algorithm  can  be  formally  given  as  follows: 

A  Updating  At  Iteration  n 

1.  Each  node  i  broadcasts  Bik  (fV£)  and  Gik  (F^)  for  each 
keO(i) 

2.  After  receiving  the  above  information  from  all  nodes,  each  node  i 
computes  the  path  lengths  and  dpP  for  each  p  ePa  such  that 

Sp°  >0  and  >  0: 


ivn 


,  ,  ^\n  <8(k  <C>  *  G(k  (F"k)} 

U,k)ep:Sp  >0  **  U 
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E 

(l,k)ep:S 


dn 

P 


>0 


I 


3. 


The  minimum  voice  and  data  distances  between  nodes 
then  computed  as  follows: 

dvn(a)  =  Min  dvn 
P£Pa  P 


and  j  are 


ddn (a)  *  Min  ddn 

PeP  a  P 

Usually  there  are  many  paths  from  node  i  to  node  j  ,  and  it  is  not 
computationally  feasible  to  find  the  length  of  every  path.  This  is  why 
in  step  2  we  require  that  only  the  lengths  of  those  paths  with  >  0 
and  Sp°  >  0  be  computed  for  voice  and  data  conversations,  respectively. 

4.  For  each  a  =  (i,j)  node  i  computes  the  voice  priority  function 
Pa(*a)  if  $vn(a)  >  0,  and  the  data  priority  function  qa ( r")  if  Sdn(a)  >  0. 


B.  Voice  Flow  Control 


5.  Let  Zj  =  dvn(a)  -  pa(£j) 


6.  % 


n+1 


(3.21) 


if  0  <  *2  *  Va  <  aij 


if  *a"  -  \  Za  -  aij 
if  *a"  -  **  Zan  1  0 


(3.22) 


where  n.  is  some  positive  scale  factor,  and  a.,  is  the  nominal  voice 

X  1 J 

packet  length  defined  earlier.  What  this  algorithm  does  is  the  following: 
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J 


(i)  if  PaUa)  >  dvn{a)  * 
(ii)  if  Pa(*a)  <  dVn(a), 


increase  a 
decrease  l 


n 

a  * 
n 

a  * 


£  <  a 
a  ij 

>  0 
a 
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where  py  is  a  voice  shortest  path.  What  this  algorithm  does  is  the 
following.  It  measures  the  deviation  e^11  of  the  length  of  each  path 
pePa  from  the  shortest  distance.  Then  it  decreases  the  flow  on  path  p 

<*  i 

by  an  amount  proportional  to  this  deviation,  and  routes  the  remaining 
voice  traffic  from  voice  conversation  acAv  along  a  shortest  path. Note 
that  Sv^n+^(a),  the  total  voice  traffic  entering  the  network  at  node  i 
and  destined  for  node  j  at  iteration  (n+1),  is  determined  by  the  voice 
flow  control  algorithm. 


E.  Data  Traffic  Routine 


Let  e 


djn  -  ddn(a) 
min  [Sdn  ,  nd  edn] 


(3.28) 


(3.29) 


12.  S 


d(n+l) 


Sd"  -Adn 
P  P 


P  *  Pa 


Sd(n+1)  (a)  -  Z*  Sd(n+1),  p  -  p*  (3.30) 

P*Pd  P 


Qn 

where  nd  is  some  positive  scale  factor,  Sp  >0,  and  pd  is  a  data 
shortest  path.  The  operation  of  this  algorithm  is  similar  to  that  of  D. 
Note  that  we  assume  that  at  each  node,  voice  packets  would  be  transmitted 
before  the  data  packets,  in  the  spirit  of  the  non-preemptive  priority  voice 
packets  have  over  data  packets.  The  convergence  properties  of  the  joint 
flow  control  and  routing  algorithm  are  described  by  the  following 


Theorem  3.2  Let  u°  be  any  feasible  point  of  (2.14)  through  (2.16),  and 
let  J(u°)  <_ JQ.  Then  for  every  positive  number  JQ,  there  exist  scale 
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factors  n^,  nr»  nv  and  nd  such  that 

Min  J(u) , 
u 

for  all  n  >_  1 

This  theorem  states  that  there  exist  n„,  n„  and  n.  such  that 

x,  r  v  a 

if  we  start  at  any  feasible  point  u°  and  apply  the  point-to-set  mapping  H 
repeatedly,  the  set  of  points  {u11}  we  obtain  converges  to  a  point  which 
minimizes  J.  The  proof  of  this  theorem  is  given  in  Appendix  B. 

3.3  Discussion 

Our  routing  algorithm  is  similar  to  the  new  ARPANET  routing 
algorithm  [26]:  both  are  shortest  path  routing  algorithms.  However, 
each  uses  a  different  metric  for  distance  computation.  In  [26]  the 
distances  are  absolute  delay.s  averaged  over  10  seconds,  but  in  our  work 
the  distances  are  the  marginal  delays  on  the  paths.  Also  in  our  work 
we  gradually  decrease  the  traffic  on  the  paths  that  are  no  longer  the 
shortest  paths  at  each  iteration,  thereby  ensuring  that  the  algorithm 
converges.  But  in  [26],  at  each  iteration  an  old  path  that  is  no  longer 
the  shortest  path  is  entirely  cancelled  and  all  traffic  is  shifted  to 
the  new  shortest  path.  This  type  of  algorithm  is  not  likely  to  converge  [2]. 
Our  routing  algorithm  is  also  similar  to  the  algorithm  proposed  by 


Lim  J(u  ) 
rW] 


where 


un  e  H(un_1)  , 


Segal!  [29].  Both  are  shortest  path  algorithms  arid  use  the  marginal 
delay  for  distance  computation.  Also  in  both  algorithms  the  computation 
is  distributed.  However,  the  two  algorithms  differ  in  a  number  of  ways. 
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First,  in  [29]  the  algorithm  is  defined  in  the  space  of  link  flows. 

Secondly,  the  update  information  in  [29]  is  propagated  upstream,  starting 
from  the  destination  node  as  follows:  Each  node,  upon  receiving  the 
update  information  from  its  downstream  neighbor,  performs  the  necessary 
computation  and  forwards  the  result  to  its  upstream  neighbors  who  in  turn 
perform  similar  operations  and  pass  the  information  further  upstream, 
until  the  source  node  is  reached.  This  mode  of  update  information 
transmission  has  the  potential  for  generating  loops  that  could  prevent 
the  upstream  nodes  from  ever  receiving  the  update  information.  Specifically, 
if  there  is  a  link  (k,£)  such  that  f^  (j)  >  0  and  f^.  (j)  >  0,  then  during 
an  update  node  k  needs  the  value  of  the  distance  X^  to  compute  X^ 
and  node  A  needs  the  value  of  Xfc.  to  compute  X^-.  Then  a  deadlock  ensues 
and  no  further  upstream  forwarding  of  the  update  information  can  take 
place.  This  necessitates  the  introduction,  as  was  done  in  [29],  of  the 
concept  of  blocking  which  ensures  that  loop-freedom  is  maintained.  In 
cur  algorithm  the  update  information  is  broadcast  and  so  we  do  not  worry 
about  loop- freedom. 

Our  routing  algorithm  is  also  similar  to  the  algorithm  proposed  by 
Gallager  [12].  However,  that  algorithm  is  defined  in  the  space  of 
fractions  of  link  flows  rather  than  path  flows  in  which  our  algorithm 
is  defined.  Also,  because  of  the  way  the  update  information  is  trans¬ 
mitted  (same  as  for  [29]), the  issue  of  blocking  is  introduced  to  maintain 
loop- freedom.  The  original  idea  of  blocking  as  a  means  of  maintaining 
loop- freedom  was  formulated  in  [12]. 

Finally,  our  flow  control  algorithm  is  similar  to  that  of 
Golestaani  [19].  This  is  because  the  two  models  are  similar,  differing 


•6 


only  in  the  fact  that  we  have  two  additional  cost  functions  for  voice 
traffic. 


A  Class  of  Priority  Functions 


We  have  defined  the  voice  and  data  priority  functions  as 

*1  * 

PmUm)  =  -  5-  V„  (O  ,  m  £  Aw 

rm  m'  g  m  m  v 


qm(rn,>  *  - 


respectively.  Golestaani  [19]  introduced  the  concept  of  priority 
functions  in  data  networks.  He  also  considered  a  class  of  priority 
functions  which  enables  us  to  prioritize  the  different  conversations 
according  to  how  severely  they  are  to  be  cut  back  when  the  network 
loading  becomes  high.  We  shaTl  consider  this  class  of  priority  functions 
in  this  section  together  with  a  modified  form  of  this  class  of  priority 
functions  for  data  conversations.  Consider  priority  functions  of  the  form: 

Vm 

m 

vm  >  0  (3.31) 

5m  Sm  >  0  (3.32 

m 

We  define  v_  and  <5_  as  the  voice  priority  index  for  voice  conversation  m, 
m  m  - c - ■* — - 

and  the  data  priority  index  for  data  conversation  respectively.  The 
quantity  am  is  the  nominal  voice  packet  length  defined  earlier.  The 
quantity  bm  is  the  nominal  input  rate  for  data  conversation  m.  Note  that 
we  do  not  know  r|j|  a  priori,  but  we  still  assign  the  nominal  rate  bm  to 
conversation  m.  Figure  3.1  is  a  plot  of  p  (t  ). 


-47- 


Figure  3.1  Voice  Priority  Function 


The  priority  indices'  relate  to  the  extent  to  which  conversations 
get  cut  back  when  the  network  loading  becomes  high.  Conversations  with 
small  priority  indices  are  cut  back  more  severely  than  those  with  higher 
priority  indices.  As  can  be  seen  in  Fig.  3.1,  packet  lengths  of  voice 
conversations  with  high  \>m  remain  fairly  constant  over  a  wide  range  of 
values  of  p  But  packet  lengths  of  voice  conversations  with  lower 

values  of  vm  vary  wildly  as  pmUm)  changes.  When  all  conversations 
have  high  priority  indices  it  is  no  longer  true  that  packet  lengths  remain 
fairly  constant.  This  is  because  when  congestion  sets  in  all  packets 
would  have  to  be  shortened  to  reduce  network  traffic.  And  since  all  the 
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conversations  have  the  same  order  of  importance,  they  would  all  then 

experience  severe  cutback  in  packet  length. 

As  we  noted  in  chapter  II,  the  nominal  voice  packet  length  am 

relates  to  the  length  of  packets  required  to  produce  good  quality  speech. 

However,  for  fixed  values  of  p  (£)  and  v  .  the  higher  the  value  of  am 

m  m  m  m 

the  higher  is  the  value  of  £  .  Since  the  amount  of  traffic  generated  by 

a  voice  packet  is  directly  proportional  to  its  length,  care  must  be  taken 

not  to  choose  such  values  of  am  that  could  indirectly  induce  congestion. 

Similar  arguments  hold  for  the  choice  of  bm- 

We  mentioned  earlier  that  we  do  not  know  rj|j  a  priori  but  still 

assign  a  nominal  rate  bm-  It  may  turn  out  that  r||J  >  bm  and  an  assigned 

rate  r,  >  b  .  Observe,  however,  that  when  r  >  bm  some  inversion  takes 
mm  mm 

place:  For  any  given  value  of  qm(rm)  <  T,  the  data  conversations  with 

lower  priority  indices  get  higher  rates  than  those  with  higher  priority 

indices.  If  this  is  the  desired  objective  of  the  network  designer,  then 

the  definition  of  qm(rm)  given  in  (3.32)  is  enough.”  However,  it  may  be 

desirable  to  allow  the  higher  priority  conversations  (i.e.  those  with 

high  priority  indices)  to  have  higher  input  rates  at  all  times.  For 

this  purpose  the  definition  of  q  (r  )  could  be  modified  as  follows: 

mm 


\(rJ  ■ 


0  <  r_  <  b_ 
m  —  m 


r  >  bm 
m  —  m 


(3.33) 


The  modified  priority  function  is  sketched  in  Fig.  3.2. 

5. 


*3  2 

\\ 


.increasing  5 


Figure  3.2  Modified  Data  Priority  Function 

For  the  class  of  priority  functions  defined  by  (3.31)  and  (3.32) 
we  have  that  for  voice  conversations  me: Av  ,  when  *  am,  Pm(Jtm)  *  1. 
When  this  happens  we  would  not  like  the  network  to  be  congested,  and 
we  would  not  like  to  unnecessarily  restrict  the  flows  to  some  arbitrarily 
small  values  since  this  is  one  way  of  avoiding  congestion.  Recall  that 
the  voice  flow  control  algorithm  operates  as  follows: 


(i) 

if 

p.  <0 

>  dvn(m)  , 

it 

increases 

C 

if  *:<  v 

(ii) 

if 

<  dvn(m). 

it 

decreases 

if  «."  >  0  ; 

m 

(iii) 

if 

p jO 

a  dvn(m). 

it 

leaves 

c 

unchanged. 

If  our  network  is  properly  designed  (and  we  have  no  reason  to  believe  it 

is  notl)  ,  it  will  operate  at  the  equilibrium  point  most  of  the  time. 

That  is,  p  (t")  =  dvn(m)  most  of  the  time.  We  would  like  to  know  what 
mm 

happens  when  PmUm)  =  1  and  the  network  operates  at  the  equilibrium  point. 
Recall  that 

-  8ik  <F™>  +  <4  <4>  • 

dvn(m)*  Min  £  d™ 

<*Pn  (l.k)cp 

Assume  that  all  links  (i,k)  ep  have  the  same  flow,  and  hence  the  same 

vn 

value  of  d  and  let  cr  be  the  number  of  links  on  the  path  peP. 

ik  m  r  r  m 

Then  we  require  that  at  the  equilibrium  point  when  Pm(^)  =  1. 

am  dik  ~  ^ •  That  1S»  dik  *  Vam.  Tiius  when  the  network  °Perates  at 
the  equilibrium  point  with  =  am,  the  voice  length  of  each  link 
(i,k)  ep  is  determined  by  the  number  of  links  on  the  path  p.  Furthermore, 
as  the  number  of  links  am  increases,  the  voice  length  d^  decreases. 

For  example,  when  =  2,  dV£  =  0.5;  and  whenam  =  5,  dV£  =  0.2.  This 
means  that  for  a  multi-hop  conversation  m,  the  condition  that  PmUm)  *  1 
at  the  equilibrium  point  implies  that  the  traffic  on  its  path  be  very 
light.  The  above  argument  for  voice  conversations  hold  for  data 
conversations.  In  sum  then,  when  we  choose  a  class  of  priority  functions 
PmUm)»  we  indirectly  define  the  nature  of  the  class  of  congestion  cost 
functions  that  are  compatible  with  it.  This  is  so  because  the  optimality 
conditions  demand  that  there  be  a  proper  scaling  between  the  congestion 
cost  functions  and  the  priority  functions. 
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CHAPTER  IV 

EXTENSIONS  TO  THE  BASIC  ALGORITHM 

4.1  Introduction 

The  functions  of  the  basic  algorithm,  the  joint  flow  control  and 
routing  algorithm,  are  the  setting  of  the  voice  packet  lengths  and  data 
input  rates  to  values  compatible  with  the  network  utilization,  and  the 
routing  of  the  voice  and  data  traffics.  In  this  chapter  we  consider  a 
group  of  protocols  for  extending  the  functions  of  the  algorithm.  In 
particular,  we  consider  a  quasi-static  protocol  that  is  used  to  reduce 
the  likelihood  of  discarding  many  voice  packets  which  may  arrive  late 
at  the  smoothing  buffers  due  to  high  network  loading.  We  also  consider 
protocols  for  dealing  with  congestion  on  a  more  dynamic  basis  than  the 
joint  flow  control  and  routing  algorithm  can  deal  with.  Finally,  we 
consider  the  place  of  the  cut-through  routing  protocol,  proposed  by 
Kermani  and  Kieinrock  [20],  in  the  integrated  network. 

Despite  the  flow  control  we  practice,  the  network  can  still  get 
congested.  This  can  be  caused  by  several  factors  including  the  following. 
The  joint  flow  control  and  routing  algorithm  is  intended  for  quasi-static 
applications  where  the  network  input  statistics  change  very  slowly. 

However,  in  practice  the  network  has  varying  input  statistics.  Between 
updates,  because  of  the  inherent  burstiness  of  data  conversations,  some 
data  conversations  may  decrease  their  desired  rates  to  values  less  than 
the  assigned  rates,  while  others  may  terminate  their  conversations  entire¬ 
ly.  Some  inactive  data  conversations,  and  voice  conversations  in  silence 
may  become  active.  If  voice  packet  lengths  and  data  input  rates  are  as¬ 
signed  on  the  basis  of  low  network  utilization,  then  the  new  active  sources 
may  generate  such  additional  traffic  that  could  drive  the  network  into  con- 


gestion.  Therefore,  the  joint  flow  control  and  routing  algorithm,  while 
reducing  the  likelihood  of  network  congestion,  is  not  a  perfect  safe¬ 
guard  against  congestion  because  it  will  not  be  able  to  track  the  chang¬ 
ing  i-nput  statistics.  To  do  this,  the  algorithm  would  have  to  be  up¬ 
dated  more  frequently,  which  requires  more  updating  protocols  and  hence 
a  reduction  in  the  effective  link  capacities  available  for  the  voice  and 
data  traffics. 

4.2  A  Direct  Method  of  Congestion  Control 

There  are  differing  views  on  what  to  do  when  the  network  gets  con¬ 
gested.  Some  researchers  propose  that  voice  packets  be  dropped  when  net¬ 
work  congestion  occurs;  and  if  congestion  persists,  then  some  data  packets 
be  dropped  next  [11].  The  reason  for  this  is  that  it  is  easier  to  drop 
voice  packets  because  the  speaker  can  easily  repeat  the  message  if  he 
does  not  get  a  response  in  time.  Moreover,  since  the  voice  packets  are 
likely  to  arrive  late  at  their  destinations  and  hence  be  discarded,  it 
is  better  to  get  rid  of  them  earlier  rather  than  waste  the  network  resources 
in  transmitting  them.  The  problem  with  dropping  voice  packets  as  a  means 
of  fighting  congestion  is  that  if  the  algorithm  that  determines  when  voice 
packets  are  to  be  dropped  is  not  dynamic  enough,  it  would  not  stop  drop¬ 
ping  packets  after  relieving  congestion. 

The  embedded  coding  scheme  of  Bially  et  al .  [3]  controls  congestion 
In  a  dynamic  manner,  and  is  adaptable  to  our  model.  Recall  that  in  that 
scheme  voice  is  encoded  into  the  more  important  packets  and  the  less 
important  packets.  The  less  important  packets  may  be  dropped  anywhere  in 
the  network  on  a  dynamic  basis,  according  to  the  network  loading  conditions. 
The  source  node  would  have  to  be  informed  of  any  packets  that  have  been 
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dropped  so  that  it  would  generate  only  the  more  important  packets.  We 
emphasize  at  this  point  that  this  scheme  is  equivalent  to  using  one  packet 
that  has  different  segments  of  different  orders  of  importance.  When  the 
network  loading  becomes  high,  the  less  important  segments  of  the  packet 
could  be  discarded  anywhere  in  the  network.  Flow  control  is  practised 
by  setting  the  voice  packet  length  to  a  value  compatible  with  the  network 
conditions  before  sending  it  into  the  network.  This  corresponds  to  drop¬ 
ping  the  less  important  segments  of  the  packet  at  the  source. 

We  would  like  to  make  a  distinction  between  flow  control  and  con¬ 
gestion  control,  even  though  the  two  terms  are  often  used  interchangeably 
in  the  literature.  For  our  purpose  we  define  flow  control  as  any  mechanism 
that  regulates  the  entry  of  traffic  Into  the  network  in  order  to  avoid 
network  overload.  In  other  words,  flow  control  is  applied  at  the  entry 
points  of  the  traffic.  Thus,  voice  packet  length  setting  and  data  input 
rate  setting  are  examples  of  flow  control.  By  congestion  control  we  mean 
any  mechanism  that  alleviates  high  network  loading  by  acting  on  the  traf¬ 
fic  already  inside  the  network.  Thus  by  this  definition,  the  embedded  cod¬ 
ing  scheme  of  Bially  et  al .  [3]  is  a  congestion  control  scheme  and  not  a 
flow  control  scheme. 

In  addition  to  effecting  flow  control  we  propose  to  practise  congestion 
control  in  the  integrated  network.  Our  approach  is  similar  to  the  embedded 
coding  scheme  [3].  After  setting  the  voice  packet  lengths,  each  voice 
packet  would  be  encoded  into  different  segments  of  different  orders  of 
importance.  As  in  [3],  when  the  network  loading  becomes  high,  the  less 
important  segments  of  the  packet  could  be  discarded  anywhere  in  the  net¬ 
work.  The  more  important  segments  reaching  the  destination  are  considered 
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capable  of  producing  an  acceptable  speech  quality,  but  the  quality  is 
lower  than  would  be  the  case  if  the  packet  was  received  with  all  segments 
in  place. 

The  difference  between  our  proposed  scheme  and  that  in  [3]  is  that 
we  first  practise  flow  control  before  practising  congestion  control.  In 
[3]  no  flow  control  is  practised.  Also  we  have  no  need  for  a  mechamism 
for  informing  the  source  node  of  the  dropped  segment  because  the  flow  con¬ 
trol  algorithm  will  take  care  of  that  by  setting  the  voice  packet  lengths 
to  compatible  values  at  the  next  iteration. 

4.3  The  Variable  Voice  Packet  Delivery  Rate  Strategy 

We  require  that  voice  packets  be  delivered  to  the  sink  from  the  smooth¬ 
ing  buffers  at  the  same  rate  8  at  which  the  voice  digitizer  emits  them  in¬ 
to  the  network.  However,  when  the  network  gets  congested,  voice  packets 
arrive  very  late  and  so  the  smoothing  buffers  would  be  delivering  packets 
to  the  sink  faster  than  new  packets  arrive  at  the  buffers.  In  this  case 
many  of  these  packets  would  be  discarded  because  of  late  arrival. 

In  order  to  reduce  the  number  of  voice  packets  that  are  discarded, 
we  would  have  to  vary  the  rate  at  which  the  smoothing  buffers  deliver 
these  packets  to  the  sink  in  response  to  the  network  conditions.  Specifical 
ly,  when  the  network  loading  becomes  high  and  packets  arrive  later  than  be¬ 
fore,  we  decrease  the  packet  delivery  to  a  value  8'.  When  the  network 
loading  becomes  light  again,  we  increase  8'  back  to  8.  In  order  to  fully 
understand  the  effect  of  changing  the  packet  delivery  rate  to  S'*  con¬ 
sider  the  delay  of  a  packet  that  enters  the  network  at  time  t  and  leaves 
the  network  (i.e.  is  delivered  to  the  sink)  at  time  d;  see  Figure  4.1. 
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x  t  +  1/6  t  +  n/B 


Figure  4.1  Additional  Delay  Caused  By  Change  In  3 

As  can  be  seen  from  the  figure,  a  packet  that  enters  the  network  at  time 
t  +  n/B  is  expected  to  be  at  the  smoothing  buffers  no  later  than  d  +  n/B', 
sustaining  a  maximum  delay  of 

(d  +  -~r)  -  (t  +  n/B  )  =  (d-t)  + 

without  being  thrown  away.  If  S  =  8',  then  this  packet  would  be  expected 
to  suffer  a  maximum  delay  of  (d-r).  But  with  3  f  6',  we  permit  the  packet 
to  suffer  an  additional  delay  of  ^  without  being  thrown  away  when 

it  arrives.  Note  that  this  additional  delay  increases  with  time.  Thus 
the  packets  that  enter  the  network  long  after  we  have  switched  the  deliv¬ 
ery  rate  from  3  to  S'  have  less  chances  of  late  arrival  than  those  that 
enter  immediately  after  the  switch.  Since  the  variation  of  voice  deliv¬ 
ery  rate  affects  the  voice  quality,  the  value  of  n(B-B’)/6B'  should  be 
such  that  Is  a  compromise  between  the  probability  of  discarding  fewer 
packets  and  the  amount  of  quality  degradation  the  sink  can  tolerate.  We 
can  define  a  priori  how  large  we  would  make  the  additional  delay 
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n(8-B')/8B' .  When  this  predetermined  value  is  reached,  we  revert  to  $ 

whether  congestion  still  persists  or  not. 

4.4  A  Congestion  Avoidance  Protocol 

The  joint  flow  control  and  routing  algorithm  uses  yn. .  as  its  esti- 

*  0 

mate  of  nl.  in  setting  the  voice  packet  lengths.  This,  as  we  have  noted, 

*  j 

is  due  to  the  fact  that  it  is  difficult  to  use  the  instantaneous  values 
of  nl.  on  a  global  basis.  This  difficulty  results  in  occasional  inci- 

•  J 

dence  of  heavy  network  loading,  which  in  turn  necessitates  such  actions 
as  the  shortening  of  voice  packets  and  the  varying  of  the  voice  packet 
delivery  rate,  as  we  recommended. 

It  is  possible  to  use  the  nj.  on  a  local  basis,  especially  in  a 
manner  that  would  enable  each  source  node  i,  which  sends  traffic  to 
destination  node  j,  to  permit  some  maximum  variation  of  the  (i,j)  voice 
traffic  from  the  expected  rate.  The  total  expected  rate  of  all  voice 
(i,j)  conversations  is  Byn.  .Jit,,  where  2*.  is  the  value  of  the  voice 

*  u  '  J  ■  J 

packet  length  determined  by  the  algorithm.  If  at  any  time  greatly 

*  J 

exceeds  youthen  the  network  may  be  tending  to  congestion  because  of 
the  increased  rate  over  and  above  the  expected  rate.  We  propose  that 
the  total  instantaneous  voice  rate  for  the  (i,j)  conversations  be  limited 
to  a  value  no  greater  than  8[yn- •  +  e(n-)]2.*-,  for  some  choice  of  e(n.  ) 

*  w  *  >J  ■  w  I J 

to  be  explained  later.  Specifically,  let  t. .(t)  be  the  instantaneous 

’  J 

value  of  5,^..  At  each  instant  that  voice  packets  are  to  be  constructed, 
note  i  makes  the  following  decision: 

(  "Tj  S 
a  ) 


,  otherwise 


(4.1) 


This  strategy  can  be  explained  as  follows.  That  we  are  dealing  with 


average  rates  means  that  we  recognize  that  fluctuations  about  the  mean 
value  exist.  But  since  voice  traffic  is  delay-sensitive,  we  cannot 
tolerate  large  upswings  in  voice  traffic.  Therefore,  we  limit  the  maxi¬ 
mum  upswing  to  +  e(n.. )]£*-.  When  traffic  exceeds  this  value,  we 

degrade  the  voice  quality  of  all  the  (i,j)  voice  conservations  through 
assigning  shorter  packets  that  would  force  the  traffic  down  to  where  we 
can  tolerate.  The  choice  of  e(n..),  the  maximum  amount  by  which  rJ., 
exceeds  before  action  is  taken,  depends  on  the  value  of  n...$pecifi- 

cally,  since  the  standard  deviation  of  nl,-  is  proportional  to  /n - . ,  we 

J  i  j 

permit  the  instantaneous  values  of  nl.  to  exceed  some  standard  deviations 

'  J 

(or  fraction  of  a  standard  deviation)  before  we  take  any  action  on  the 
traffic.  Thus  functions  of  the  form 


e  (n^j)  =  k  /no¬ 

where  k  >  0  is  a  constant,  are  appropriate  for  this  purpose. 

4.5  A  Strategy  For  Admitting  New  Conversations 

As  the  joint  flow  control  and  routing  algorithm  shows,  we  make  small 
routing  changes  at  each  iteration  in  order  to  ensure  that  the  algorithm 
converges.  Therefore,  anv  decision  to  admit  new  conversations  into  the 
network  should  be  made  in  such  a  way  that  large  changes  are  not  made  in 
the  flow  pattern.  For  all  conversations  in  the  network  let 

=  the  minimum  acceptable  voice  packet  length 
rQ  a  the  minimum  acceptable  data  input  rate 
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That  is,  if  a  voice  conversation  is  assigned  a  packet  length  less  than 

it  opts  out  of  the  network  rather  than  transmit  at  the  assigned  rate. 

A  similar  situation  holds  for  r .  For  each  voice  conversation  a,  let 

o 

X  =  some  threshold  packet  length  (whose  meaning  we  shall  explain 

a 

shortly) 

d^(a)  =  some  threshold  voice  path  length 

d^n(a)  =  the  path  length  of  the  longest  path  that  voice  conversations 
P* 

a  are  using  at  iteration  n. 

Similarly,  for  each  data  conversation  a,  let 
ra  3  some  threshold  input  rate 

a 

d^(a)  *  some  threshold  data  path  length 

d^n(a)  *  the  path  length  of  the  longest  path  that  data  conversations 
P* 

a  are  using  at  iteration  n. 

Then  the  admission  policy  we  propose  is  the  following: 

1.  If  d][n(a)  <  d]((a)  and  Z1?  >  jl  at  iteration  n,  then  and  only  then 

P*  —  o  a  —  a 

will  node  i  admit  one  new  voice  (i,j)  conversation,  where  a  *  (i,j), 
during  period  n. 

2.  If  d^n(a)  <  d^(a)  and  r"  >  r  at  iteration  n,  then  and  only  then  will 

P*  —  o  a  —  a 

node  i  admit  one  new  data  a  conversation  during  period  n. 

This  admission  policy  can  be  explained  as  follows.  We  would  not 
admit  any  new  voice  (i,j)  conversation  when  the  old  (i,j)  conversations  (i.e. 
those  already  in  the  network)  have  a  delay  greater  than  some  threshold 
d^(a).  Also  we  would  not  admit  any  new  voice  ( i , j )  conversations  when  the 
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old  conversations  have  packet  lengths  that  are  very  close  to  the  minimum 
acceptable  value  Therefore,  we  must  ensure  that  both  the  delay  per¬ 
formance  and  quality  of  all  old  a  conversations  are  of  accepable  standards 
(which  can  be  set  a  priori)  before  we  admit  any  new  a  conversation.  And 
the  reason  for  restricting  the  admission  to  at  most  one  voice  conversation 
per  source-destination  pair  per  period  is  to  make  small  changes  to  voice 
traffic  at  a  time  as  much  as  possible.  The  same  arguments  hold  for  the 
data  conversations. 

An  indication  of  how  to  choose  can  be  seen  from  Figure  4.2.  As 

the  figure  shows,  J.  is  close  to  the  nominal  packet  length  a... 

a  i  j 


h 

l 


o 


1/2  (a..  +  £0) 


a.. 


Figure  4.2  Choice  of  Threshold  Packet  Length  JTa 


Not"  that  the  newly  admitted  voice  conversation  would  be  assigned  the 
current  value  of  packet  length  and  routed  along  the  shortest  path.  Similarly, 
for  data  conversations,  the  new  a  conversation  would  be  assigned  the  cur¬ 
rent  input  rate  of  the  a  conversations.  If  the  rates  are  translated  in¬ 
to  window  sizes  (see  section  4.7),  then  the  new  data  conversation  would 
be  assigned  the  current  window  size,  and  routed  along  the  shorest  path. 

4.6  Constrained  Voice  Packet  Lengths 

As  mentioned  in  the  previous  section,  assume  we  have  a  minimum 
acceptable  voice  packet  length  lQ.  Furthermore,  let  £Q  belong  to  a  set 
L  of  admissible  voice  packet  lengths: 
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L  =  (0,  2>0>  A-j  3jj}  (4.2) 

where  0  <  <  tj  <  ...  <  a.^.  Let  the  "raw  packet  length"  obtained  from 

the  flow  control  algorithm  be  2.  :  and  let  l  satisfy  the  condition 
\-l  <  ^a  <  where  ^  £  L.  Then  we  are  required  to  set  l  to 

/v 

some  value  £  e  L. 

a 

The  problem  with  this  type  of  extra  constraint  is  that  it  could  up¬ 
set  the  convergence  of  our  algorithm.  One  way  to  get  around  this  problem 
is  to  divide  the  voice  (i,j)  conversations  arbitrarily  into  two  groups,  A  and 
B.  Then  we  use  the  following  assignment  rule: 

for  group  A 

&k_1  for  group  B  (4.3) 

where  the  dividing  line  for  group  membership  should  be  drawn  in  such  a 
way  as  to  make  the  net  traffic  remain  close  to  the  level  given  by  the 
flow  control  algorithm.  That  is, 

ynaAa  *  |A|*k  +  |B|i|c_1  (4.4) 

where  |X|  is  the  cardinality  of  set  X.  In  this  way  we  avoid  the  over- 

A 

load  caused  by  letting  all  conversations  be  assigned  2a  =  2k,  and  the 

A 

under-utilization  caused  by  assigning  i,g  *  2-k_^  to  all  conversations. 
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4.7  Implementation  of  Data  Flow  Control 

We  have  discussed  how  flow  control  can  be  effected  in  the  integrated 
voice  and  data  network:  Voice  packet  lengths  and  data  input  rates  are  set 
in  response  to  the  network  conditions.  For  voice  conversations  we  assume 
that  we  can  find  a  voice  digitizer  which  generates  packets  of  appropriate 
lengths  at  a  constant  rate  of  &  per  second  whenever  a  speaker  is  in  talk- 
spurt  [3].  Thus  if  such  a  hardware  is  found,  voice  flow  control  can  be 
easily  implemented.  For  data  conversations  the  input  rates  are  controlled 
and  the  packet  lengths  are  assumed  to  be  independent  of  network  control. 

Assume  the  expected  data  packet  length  is  r  bits  and  the  optimal  input  rate 

★  ★ 
of  data  conversation  m  is  r.  Then  in  order  to  implement  this  rate  of  r 

m  r  m 

bits/sec.,  the  packets  of  this  conversation  would  be  sent  into  the  network 

★ 

at  an  average  rate  of  rm/r  per  second.  As  we  mentioned  earlier,  that  we 

have  effected  flow  control  does  not  guarantee  a  congestion- free  operation. 

Often  statistical  fluctuations  arise  that  cause  the  flow-controlled  network 

to  be  congested.  Thus  a  strict  adherence  to  sending  packets  into  the  network 

* 

at  an  average  rate  of  rm  /r  may  not  be  a  good  strategy,  especially  when 

congestion  sets  in.  What  is  required  then  is  a  scheme  that  will  not  only 

★ 

restrict  the  input  rates  to  the  optimal  values  rm  but  also  respond  dynami¬ 
cally  to  the  network  conditions  between  updates. 

A  common  approach  to  effecting  flow  control  in  a  data  network  is 
the  so-called  window  strategy  [6,15].  In  this  strategy  each  source  node 
i  keeps  the  number  of  unacknowledged  (or*  outstanding)  packets  in  the 
network,  of  each  conversation  destined  for  node  j,  below  a  given 
number  w. .  called  the  window  size  (or  window  width).  Thus  a  new  packet 

*  J 

cannot  be  sent  into  the  network  if  the  number  of  outstanding  packets  of 
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conversation  (i,j)  at  that  time  is  equal  to  w...  This  is  a  dynamic  scheme 

*  J 

because  when  the  network  is  heavily  loaded,  packet  acknowledgements  (ACKs) 

arrive  late  and  hence  the  frequency  of  packet  entry  into  the  network 

decreases.  On  the  other  hand,  when  the  network  loading  is  light,  ACKs 

arrive  faster  and  the  frequency  of  packet  entry  into  the  network  increases. 

One  question  we  may  ask  at  this  point  is:  Given  an  optimal  input 
★ 

rate  r  ,  can  we  set  a  window  size  wm  such  that,  by  ensuring  that  the 

number  of  outstanding  data  packets  in  the  network  is  no  greater  than  wm, 

★ 

we  force  the  data  conversation  tn  to  transmit  at  rate  r  ?  It  turns  out 

m 

that  under  certain  conditions,  which  we  shall  give  shortly,  the  answer 
to  this  question  is  yes.  Golestaani  [19]  shows  that  if  for  each  active 
conversation  (i,j)  (i.e.  one  that  is  not  silent)  there  are  packets  waiting 
to  enter  the  network  so  that  as  soon  as  a  new  ACK  arrives  at  node  i  a 
new  packet  enters  the  network,  there  is  a  unique  correspondence  between 
the  set  of  outstanding  packets  in  the  network  and  the  set  of  input  rates 
of  these  conversations.  The  implication  of  this  statement  is  the  following. 
Since  the  window  strategy  is  a  dynamic  scheme,  we  can  practise  it  as  an 
indirect  enforcement  of  the  flow  control  algorithm  for  those  conversations 
that  meet  the  necessary  conditions. 

The  motivations  for  practising  the  window  strategy  in  a  data  net¬ 
work  carry  over  to  the  integrated  network.  However,  the  window  strategy 
cannot  be  applied  to  voice  conversations.  This  is  not  only  because  we 
set  the  voice  packets  to  the  desired  lengths  and  strictly  send  them  into 
the  network  at  a  constant  rate  of  3  per  second  but  also  because  the 
assigned  window  size  may  be  less  than  the  number  of  voice  packets  generated 
during  a  talkspurt.  This  would  lead  to  loss  of  continuity  in  the  conver  - 
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sation,  an  undesirable  feature  in  any  voice  conversation. 

In  deriving  the  relationship  between  r^  and  w^  Golestaani  [19] 
defined  a  function  t^F^),  which  is  the  average  delay  per  packet  on 
link  l  including  processing  and  queueing  delays  at  the  input  to  link  i. 
He  assumed  that  t^F^)  is  a  convex  increasing  and  differentiable  function 
of  F^,  the  aggregate  flow  on  link  £.  In  our  model  the  link  delay  of  a 
data  packet  is  not  only  dependent  on  the  aggregate  link  flow  but  also 
on  what  proportion  of  that  link  flow  is  voice  traffic.  For  example, 
consider  a  link  flow  of  one  unit.  If  the  voice  traffic  is  0.3  units, 
then  the  data  packet  delay  would  be  smaller  than  the  case  where  the  voice 
traffic  is  0.7  units.  In  the  latter  case  a  data  packet  is  likely  to  wait 
a  long  time  as  newly  arriving  voice  packets  receive  service  before  it. 

Thus  in  our  analysis  of  the  window  strategy  for  the  integrated  network 
we  define  t£(FA,  F^)  to  be  an  increasing  and  differentiable  function  of 
both  F^  and  F^  ,  the  voice  traffic  on  link  l.  The  effect  of  the  presence 
of  voice  packets  on  the  data  packet  delay  is  to  modulate  the  link  delay 
according  to  what  proportion  of  the  link  traffic  is  voice  traffic.  The 
complete  relationship  between  rm  and  wm  can  be  obtained  as  follows.  Let 

r  =  the  input  rate  of  data  conversation  m 
m 

r  =  the  expected  data  packet  length 

w„  *  the  window  size  of  conversation  m 
m 

am  =  the  average  delay  of  ACKs  for  conversation  m.  am  is  assumed 
to  be  independent  of  the  network  conditions  if  the  ACK 
traffic  has  priority  over  other  traffics. 
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q^  =  the  fraction  of  conversation  m  on  link  l.  For  unifilar 

flow,  such  as  in  virtual  line-switched  systems,  q  „  is 

mil 

either  1  or  0,  depending  on  whether  2,eRd  or  not. 

m 

Then  from  Little's  formula  [24]  the  number  of  packets  of  conversation  m 
outstanding  in  the  network  at  any  time  is  given  by 


(4.4) 


★  d 

For  those  conversations  with  rm  <  rm  ,  (4.4)  offers  a  nice  way  to  set  the 
window  size.  By  setting  wm  to 


w 


m 


V  VFr 


(4-5) 


we  force  conversation  m  to  buffer  the  remaining  data  so  that  whenever 

an  ACK  arrives  there  is  a  packet  waiting  to  be  sent  into  the  network. 

*  d 

However,  for  those  conversations  with  r  >  r  ,  the  allocated  window  size 

m  m 

wm  >  um*’  1S*  the  allocated  widow  size  is  greater  than  the  number  of 

outstanding  packets  in  the  network  at  any  time.  For  such  conversations, 

arriving  packets  would  tend  to  enter  the  network  earlier  than  those  of  the 

previous  category  because  there  is  a  great  likelihood  that  at  any  time 

the  number  of  packets  outstanding  is  less  than  the  window  size. 

To  summarize,  the  window  strategy  enables  us  to  partition  the  data 

*  d 

conversations  into  two  groups:  those  with  r  <  r  and  for  whom  wm  =  um, 

3  r  mm  mm 

*  d 

and  those  with  r  >  r  and  for  whom  wm  >  u  In  the  former  group  packets 
are  likely  to  be  buffered  in  wait  for  ACK  arrival  because  the  window  is 
likely  to  be  closed  (i.e.  the  number  of  outstanding  packets  equals  the 


window  size)  when  they  arrive.  In  the  latter  group  packets  are  likely 
to  enter  the  network  as  they  arrive  because  the  window  is  likely  to  be 
open  when  they  drrive. 

If  r*  <  rd  for  all  m,  we  show  below  that  when  the  window  size  is 
m  m 

★ 

set  by  (4.4)  there  is  a  unique  correspondence  between  rm  and  wm  [19J. 
First,  we  note  that 


a  £V  *  Ft 


Therefore 


3 t 

am  +  22  V1*  +rm  £  q^aFT*  k  =  m 


,3wm  _ 


rm  ^  qk*  W;  ’ 


k  f  m 


Let  T  be  a  diagonal  matrix  with  elements  given  by 

Tmm  *  am  +  ^  V  h  (fV  'V 

Q  =  a  matrix  with  elements  q^ 

R  =  a  diagonal  matrix  with  entries  =  rm 

ni 

M  *  a  diagonal  matrix  with  entries  M  *  x-f- 

*  l 


a  matrix  with  entries  Pmk  =  r  p- 
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Then  (4.6)  can  be  written  as 

P  =  T  +  RQMQT  (4.7) 

Post-multiplying  by  R  we  obtain 

PR  =  TR  +  RQMQTRT  (since  RT  =  R) 


If  we  assume  that  rm  >  0  for  all  m,  then  since  T,  R,  and  M  are  positive 
diagonal  matrices,  PR  is  a  positive  definite  matrix,  and  hence  is  invertible. 
Since  (PR) “ ^  *  R"^  P“\  and  since  R"1  exists,  then  P"^  exists.  This 

means  that  incremental  changes  in  window  sizes  lead  to  uniquely  defined 
changes  in  data  input  rates. 

In  the  way  window  strategies  are  generally  used  in  practice,  the 

window  size  of  each  conversation  remains  fixed  throughout  the  dura ton  of 

the  conversation.  However,  the  unique  correspondence  between  r  and  w 

mm 

suggests  that  wm  could  be  adjusted  in  response  to  the  network  conditions. 

3  J  *  gj* 

Specifically,  when  >  0  we  need  to  decrease  w^;  and  when  gp— <  0  we 

'  m  m 

need  to  increase  wm  [13].  One  problem  is  how  to  effect  this  window  size 

adjustment.  A  possible  method  is  to  change  wm  by  an  amount  corresponding 

to  the  desired  change  in  rm  [14];  namely 


“m  +  tt  (F  Fj)  (4.8) 

<  1 

t 

The  problem  with  this  approach  can  be  seen  in  the  following  example. 
Consider  a  simple  network  consisting  of  one  link  and  two  nodes  and  n£. 
+  This  example  is  due  to  Professor  Gallager 


Ar, 


A  w, 


m 


m 


Two  sources  transmit  at  rates  r-j  and  r 2  from  nocle  ni  as  shown  in 
Figure  4.3.  If  we  assume  that  Fv  =  0,  then 


«i  *  r  «ri *  r2' 

w2  -  ^  t(r,  +  r2) 


Figure  4.3 


(0 


To  first  order 


Aw-j  = 

Ari  rit 

_L  t  +  — 
r  r 

[Ar1  +  Ar2] 

(ii) 

Aw2  = 

Ar?  r2t 

— —  t  +  — 

r  r 

Ca»*i  +  &r2] 

(iii) 

that  r7 

=  0,  r2  >  0;  and 

¥-  ■  -i* 

3r7 

ill  = 

3r2 

-2 

so  that  Aw]  -  j  Aw2*  Then  from  (ii)  and  (iii)  we  have  that 


Aw, 


Ar-jt 


r„t 


AWr 


Ar~  '  ’2‘  A 

t-  r*  j.  ^  i  +  — —  Ar, 


=  -A-  [t  +  r,t  ]  + 


(iv) 

(v) 
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From  (iv) 


a  aw  ^ 

~5i r 


since  Aw^  =  1  ^w2  •  T^e>'e^ore»  ^rom  (v)  we  obtain 


Aw2(2t  -  rgt  )  = 


2t(t  +  r2t  ) 


2t(t  +  r„t  ) 

Aw,  =  - r~  Ar, 

*  T(2t  -  r2t  )  * 

t  \ 

If  we  assume  r2t  »  t,  then  r2t  >  2t  and 
Aw2  =  -k  t  Ar2  ,  when  k  >  0 

That  is,  as  we  increase  the  window  size  we  decrease  the  rate.  Thus  this 
method  is  good  for  adjusting  the  window  size  of  inactive  conversations 
but  not  so  for  active  conversations,  when  it  is  applied  simultaneously  to 
the  two  groups.  Increase  in  the  window  size  of  the  former  is  accompanied 
by  an  increase  in  the  input  rate;  but  in  the  latter  it  is  accompanied  by 
a  decrease  in  the  rate. 


Consider  the  following  approach.  Let 


A  30* 

Awm  =  -nrm5F- 
m 


(4.9) 


where  n  is  some  positive  scale  factor.  As  pointed  out  in  [13],  this 

method  of  window  size  adjustment  has  some  nice  properties.  When  the 

★ 

window  size  changes  by  Awm  then  to  first  order  the  change  in  J  is 
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AW 


m 


"n?9wm  ^  3rm  *  fr0m 

mm  m 


m  3w  m  3w  3r 
111  m  mm 

-  f  (Vw  .  o*)  [TR  +  RQMQTRJ  (7w  •  J*)T 


f  (Vw  •  J*)  PR  (7W  •  J*)T 


(4.10) 


*  * 

where  Vw  •  J  is  the  gradient  of  J  with  respect  to  w.  Since  PR  is 

★  * 

positive  definite,  AJ  <  0  if  J  f  0.  Thus  this  method  leads  to  a 
descent  algorithm. 

Note,  however,  that  this  method  cannot  work  for  inactive  conversations 

(whose  rm  =  0).  And  since  the  condition  for  PR  to  be  positive  definite  is 

*  d 

that  all  conversations  have  wm  =  u  ,  or  equivalently  rm  <  r  ,  this 

mm’  ^  •’mm 

*  j 

method  is  also  not  applicable  to  conversations  with  r  >  r“.  We  can 
then  make  window  size  adjustments  in  the  following  manner:  For  inactive 
conversations  which  may  demand  a  larger  window  size  when  they  become  active. 


where  A^m  is  their  desired  change  In  rate.  And  for  active  sources  whose 
★  d 

allocated  rates  rm  <  r“  ,  use 


Aw 


m 


-  nr, 


m 


3r, 


m 
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Generally  the  window  size  obtained  from  (4.5)  is  not  an  integer. 

In  practice,  however,  fractional  window  sizes  cannot  be  easily  implemented. 

A 

One  task  then  is  to  construct  integer  window  sizes  wm  given  raw  values 
wm>  It  has  been  suggested  [19]  that  the  parameter  am  be  adjusted  in 
order  that  the  value  of  w^  in  (4.5)  is  an  integer.  Thus  if  an  ACK 

i 

encounters  a  delay  a  ,  it  is  further  delayed  by  an  amount  am  -  am> 

I 

where  am  is  the  value  of  that  makes  (4.5)  an  integer.  The  manipulations 
of  the  ACK  delay  in  the  manner  we  have  described  above  to  achieve  integer 
window  sizes  means  rounding  wm  off  to  the  next  higher  integer.  This  in 
turn  means  permitting  more  traffic  into  the  network.  This  is  thus  likely  to 
upset  the  convergence  of  our  flow  control  algorithm.  To  avert  this 
potential  problem  we  make  the  following  proposal.  For  each  source-desti¬ 
nation  pair  (i,j),  node  i  will  partition  the  data  conversations  into 
two  groups,  A  and  B.  Members  of  group  A  will  be  the  most  important  data 
conversations  while  members  of  group  B  will  be  the  less  important 
conversations.  Members  of  group  A  will  get  the  window  size  allocation 

A  A 

wm  =  rwm1 >  and  members  of  group  B  will  get  the  window  size  wm  *  |_wmJ; 
where  the  partition  is  performed  in  such  a  way  that 

in.,  w  =  | A |  TwJ  +  |B|  Lw_J 
ij  m  11  m  1  1  mJ 

where  Txl  =  the  smallest  integer  greater  than  or  equal  to  x; 

LxJ  =  the  largest  integer  less  than  or  equal  to  x. 

The  exact  demarcation  line  for  group  membership  would  depend  on  where  the 
above  approximate  equality  is  best  achieved. 


4.8  Cut-Through  Routing  of  Voice  Conversations 


We  have  emphasized  the  need  to  route  the  voice  conversations  fast 
enough  for  them  to  be  available  at  the  destination  node  when  they  are 
required  for  delivery  to  the  sinks.  We  have  also  noted  that  voice 
conversations  are  not  particularly  sensitive  to  errors.  These  two 
factors  combine  to  make  the  cut  through  routing  [20]  a  good  strategy  for 
the  voice  conversations.  The  cut-through  routing  protocol  was  originally 
defined  for  data  networks  and  works  as  follows.  When  a  packet  arrives 
at  node  i  and  an  output  link  (i,k)  is  selected  for  it,  then  after 
receiving  the  packet  header  the  packet  transmission  starts  immediately 
without  link  (i,k)  waiting  to  receive  the  entire  packet  and,  therefore, 
without  checking  for  errors,  if  link  { i , k)  is  free.  Only  if  link  (i,k) 
is  busy  would  it  be  necessary  to  buffer  an  arriving  packet,  as  in  normal 
packet  switching  operation.  In  this  way  the  delay  associated  with  buffering 
to  receive  the  entire  packet  and  check  for  errors  in  front  of  an  idle  link 
is  eliminated.  A  cut  is  said  to  have  been  made  if  the  packet  successfully 
goes  through  the  link  without  waiting  to  receive  the  packet  before  trans¬ 
mission  starts. 

The  aim  of  this  section  is  to  point  out  the  possible  gains  that 
could  be  obtained  from  applying  the  cut-through  routing  protocol  to  voice 
conversations  in  the  integrated  network.  We  will  not  go  into  detailed 
analysis  of  the  scheme;  the  interested  reader  is  referred  to  [20].  We 
will  merely  make  a  somewhat  qualitative  comparison  between  the  cut- through 
routing  and  the  normal  packet-switched  operation.  We  distinguish  between 
two  types  of  cuts:  a  perfect  cut  and  a  partial  cut.  A  perfect  cut  is  the 
type  we  described  above;  where  a  packet  arrives  at  a  link  which  is  idle  and 
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"cuts  through1'  without  having  to  wait.  When  a  packet  arrives  at  an  output 
link  that  is  busy,  it  waits  until  current  service  is  completed  whence  it 
receives  service  if  there  is  no  other  packet  in  queue.  If  the  remaining 
time  to  completion  of  the  service  of  the  current  packet  is  less  than  the 
time  to  receive  the  entire  packet  and  there  is  no  other  packet  in  queue, 
then  the  transmission  of  the  packet  can  start  while  part  of  the  packet 
is  still  being  received.  This  type  of  cut-through  is  called  a  partial 
cut  since  some  waiting  is  involved,  but  the  waiting  time  is  less  than  the 
time  to  receive  the  entire  packet.  This  virtual  waiting  time  may  also  be 
equal  to  the  time  to  receive  the  entire  packet  if  the  current  service 
ends  just  as  the  last  bit  of  the  packet  is  received.  In  [20]  only  perfect 
cut-through  routing  is  considered. 

We  present  an  approximate  analysis  of  the  cut- through  routing  to 
show  the  type  of  saving  in  packet  delay  that  can  be  derived  from  practising 
the  routing  strategy.  Assume  there  is  negligible  propagation  delay  on 
each  link  and  let 

s^  *  the  total  service  time  of  a  tagged  packet  on  link  (i,k) 

Then  s.^  is  made  up  of  two  parts: 

sik  “  wik  +  Tik 

where  =  the  waiting  time  of  the  packet  on  link  (i,k) 

=  the  transmission  time  of  the  packet  on  link  (i,k) 

Without  cut- through  routing,  each  packet  encounters  both  a  waiting  time 
and  a  transmission  time  on  each  link:  Even  when  no  other  packet  is  in 
service,  the  tagged  packet  must  be  completely  received  before  transmission 
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can  start.  Thus  the  total  system  time  of  the  tagged  packet  when  cut- 
through  routing  is  not  practised  is 


2  Cw.k  +  T  ..  ) 
(i,k)cp  lk  lk 


(4.11) 


where  p  is  the  path  of  the  tagged  packet.  With  cut-through  routing, 
for  links  (i,k),  (k,t) e  p,  the  waiting  time  at  link  (k,£)  starts  when 
transmission  starts  at  link  (i,k),  Thus  except  at  the  output  link  of  the 
destination  node  where  the  entire  packet  may  possibly  be  completely  received 
before  being  delivered  to  the  sink,  the  transmission  time  at  any  inter¬ 
mediate  link  is  absorbed  by  the  waiting  time  at  the  successor  link  on 
the  path.  Therefore,  the  system  time  of  the  tagged  packet  is 


s„  *  2  w.k  +  t.  (4.12) 

P  0,k)cp  lk  f 

where  xf  -  the  transmission  time  at  the  final  link.  From  (4.11)  and 
(4.12)  we  see  that  the  system  time  is  reduced  by  the  transmission  times 
on  all  but  the  final  output  link.  Note  that  in  practice  the  values  of 
the  w^  used  in  (4.11)  may  not  be  equal  to  those  used  in  (4.12)  because 
of  the  fact  that  perfect  cuts  might  be  made  in  some  links.  Under  such 
conditions  the  reduction  in  delay  would  be  greater  than  we  stated  above. 

The  tradeoff  in  the  cut-through  routing  strategy  is  between  delay 
and  errors.  A  packet  that  makes  a  cut  at  each  link  could  arrive  at  the 
destination  with  such  a  serious  error  that  would  necessitate  its  being 
discarded.  In  this  situation,  we  would  have  wasted  the  network  resources 
in  transmitting  that  packet.  Without  cut-through  routing,  such  a  packet 
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would  have  been  discarded  earlier,  if  the  error  proved  to  be  uncorrect- 
able,  without  the  network  resources  being  wasted  on  the  packet.  However, 
there  are  certain  errors  that  are  tolerable  in  voice  communication.  This 
means  that  we  can  practice  cut-through  routing  of  voice  traffic  in  the 
integrated  network.  We  hope  that  most  of  the  errors  we  would  encounter 
in  the  network  would  be  the  tolerable  type.  Under  this  condition  then 
the  cut- through  routing  scheme  would  provide  some  assistance  in  making 
the  voice  packets  meet  the  stringent  delay  requirements  of  the  integrated 
network. 
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CHAPTER  V 

A  CLASS  OF  INTEGRATED  NETWORK  MODELS 


Our  network  model  belongs  to  a  broader  class  of  models  for  integrating 
n  different  types  of  traffics  that  have  different  degrees  of  delay 
sensitivity  onto  one  network.  In  this  type  of  network  problem  it  is  very 
necessary  to  ensure  that  the  network  is  not  heavily  loaded  with  traffic 
from  the  very  delay-sensitive  traffic  type.  Assume  we  index  the  traffic  types 
in  a  decreasing  order  of  delay  sensitivity.  That  is,  traffic  type  1 
is  the  most  delay-sensitive,  followed  by  traffic  type  2,  and  so  on  with 
traffic  type  n  the  least  del  ay- sensitive.  Then  we  assign  the  highest 
(non-preemptive)  priority  to  traffic  type  1,  followed  by  traffic  type  2, 
etc.  The  reasons  for  practising  non-preemptive  priority  are  the  same  as 
those  advanced  in  chapter  I;  namely,  to  avoid  the  confusion  that  could 
arise  if  the  part  of  an  item  that  had  been  serviced  before  a  preemption 
occurs  is  interpreted  as  network  errors,  and  to  avoid  using  the  network 
resources  to  transmit  something  that  would  later  be  discarded.  Let 


the  expected  total  traffic  on  link(i,k)  of  traffic 
type  j,  j  =  1 ,...,n. 


+  F 


3 

ik 
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F1k  +  F?k  *■  Fix  +  F?k 

n-l  . 

£  Fik 

j=l 


\  ^  =  the  aggregate  traffic  on  link  (i,k) 

j=1 

We  define  the  cost  functions  (F^)  and  Gjk  (Fik)  as  follows: 

Bik  ^Fik^  =  the  con9est-i°n  cost  t0  link  (i»k)  of  limiting 
the  traffic  from  traffic  types  1  through  j  to 

Fik  ’  J  *  1  »2 . n-l . 

Gik  ^Fik^  *  the  congestion  cost  to  link  (i,k)  of  limiting 
the  aggregate  traffic  to  Fik. 

For  ease  of  analysis  we  assume  that  these  cost  functions  are  convex 
increasing  and  twice  differentiable  functions  of  their  respective  arguments. 
For  the  purpose  of  ensuring  that  the  network  is  not  heavily  loaded  with 

the  very  delay-sensitive  traffics,  we  require  that  8^  cuts  off  earlier 

2  3 

than  B^»  which  in  turn  cuts  off  earlier  than  B.^,  and  so  on  with  Gik 

cutting  off  the  last  but  before  the  link  capacity  C^.  By  defining 

appropriate  rate  limitation  functions  as  we  did  in  Chapter  II,  we  ensure 

that  no  traffic  type  is  completely  excluded  from  using  the  network. 

For  the  purpose  of  illustrating  how  the  congestion  cost  functions  cut  off, 

we  present  typical  plots  of  these  functions  on  the  same  set  of  axes  in 

Figure  5.1.  In  the  plot,  traffic  types  1  and  2  are  the  very  delay- 


sensitive  types,  while  types  3  through  n  are  less  delay-sensitive. 


Figure  5.1  Congestion  Costs  of  Link  Traffics 


A  more  general  communication  network  application  of  this  genera¬ 
lized  model  than  the  type  we  consider  in  this  thesis  is  the  following. 
We  can  partition  all  communication  network  traffic  into  three  different 
classes  [9]: 

Class  I  :  characterized  by  long  messages  that  require  continuous 

real-time  delivery;  e.g.  voice  facsimile,  video. 

Class  II:  characterized  by  short  discrete  messages  that  require 

near  real-time  delivery;  e.g.  interactive  data. 

Class  III:  characterized  by  long  messages  that  require  neither 
continuity  nor  immediate  delivery;  e.g.  bulk  data. 

To  integrate  these  three  traffic  types  onto  one  network,  we  have 
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to  recognize  the  need  for  classes  1  and  II  traffics  to  suffer  less 
delay  than  class  III  traffic.  This  necessitates  assigning  the  highest 
(non-preemptive)  priority  to  class  I  traffic,  followed  by  class  II  traffic, 
and  class  III  traffic  is  the  low  priority  traffic.  Let 


ik 

F?k 


F?k 


Ffk 


p|k  =  the  expected  total  class  I  traffic  on  link  (i>k) 

the  expected  total  class  II  traffic  on  link  (i,k) 

the  expected  total  class  III  traffic  on  link  (i,k) 

F1  +  F2 
ik  hik 

C1  .  c2  4.  c3 

Fik  +  Fik  +  Fik 


Then  applying  our  argument  on  the  nature  of  the  congestion  cost  functions 
we  define  Bjk  (Fjk),  B?k(F?k)  and  Gik  (Fik)  as  shown  in  Figure  5.2- 
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APPENDIX  A 


In  this  appendix  we  give  the  proof  of  theorem  3.1  which  deals  with 
tne  necessary  and  sufficient  conditions  for  optimality  of  the  solution  co 
(2.13)  subject  to  (2.14)  through  (2.16).  We  prove  the  sufficiency  part 
of  the  theorem  here;  the  proof  of  the  necessity  part  of  the  theorem  is 
given  in  Appendix  B. 

Lemma  A.l  For  any  feasible  point  P*  =  {fj*  (j),  f^Cj).  £*.*>  r*-}  which 
satisfies  the  conditions  of  theorem  3.1,  and  any  other  feasible  point 
P  =  (f^k(j),  f^k(j),  r^.}  the  following  inequalities  hold: 


(i) 


i,j,k 


> 


Yn.. . 
U 


l  x 

ij  ij 


(id  l«;„  (F*fk)  ^(j) 

i ,  j ,  k 


i  r 

i  J 


mij  rij  M i j 


with  equality  if  P 


★ 

P 


Proof  From  (3.1)  we  have  that 
-v* 


Bik<Fik)  +  Gik'F1k>+  XkJ  2-  X1j 


(*)  with  equality 


if  f^k(j)  >  0.  Multiply  both  sides  of  (*)  by  f ^ k ( j )  and  sum  over  i , j , k : 

£ce;k  *  s;k  (f;k)]f;k(3)  *  £  *kj  f?k  u>  >  £  xu  4u> 

i,j,k  i.j.k  i,j,k 
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with  equality  if  f^(j)  =  (j).  Thus 


.E  cbJ 

i »  j  » k 


ik 


t^it) +  s;,  (f’j]  (j)  > 


’ik  v'ik 


ik 


1  >  J  >  K 


'ij  fIk<J> 


E 

i  ij  »k 


kj 


f?k«> 


i.j  J 

i?fik^  "  k5 

«■ 

II 

cEf-k(j)  -  E  fjf  an 

k  5. 

*  2  xii 

8yn..  *.. 

ij  i  j 

where  the  last  equality  follows  from  (2.14).  Similarly,  from  (3.2)  we  can 
show  that  inequality  (ii)  holds,  thus  completing  the  proof  of  the  lemma. 

Now  let  there  be  and  satisfying  (3.1)  through  (3.4)  at 
a  feasible  point  P*  =  ff-k(j),  fd*kCj) .  t*j.  r*  )  .  For  any  other  point 
P-  ^k(j),  5-^.  r^}  and  a,  where  0  <_  a  £  1 ,  define 

fik (J)  •  (1  ■  a>  f(k  <■*> +  °  fik(j) 
ffk  (j)  ■  0  -  a)  ffk  (j)  +  a  fjk(j) 

(1  -  a)  tjj  *  at,j 
(1  -  a)  r*^  *  ortJ 

Then  0(a)  ■  E  [Bfk(F™)  ♦  S,k<F?k)]  *  E  [yn^  vfj  («?,)  ♦  m(J 


is  a  convex  function  of  a  .  Therefore, 


0(1) 


■  J(0)  if51 


a  =  0 


X  cBik  +  <4  C4CJ1  - 


i  .3  ,k 


+  H  Yn,.  v!.U*.)  [*•  •  -  1*2 


i  .3 


tj  iJ  iJ  u  ij  tj- 


+  E.  4  <4>  [f?k°>  •  4  (J)] 


i  »j  ,k 


ik 


*  £  »14  4  4j>  Cru  -  r*j] 
i  X  xu  -  »*j]  *  2>  "u  p(j(<j)  r* 

1  > J  i  ,3 

+  £  “u  "ij  [rij  -  4] +  £  mi3  <u<rij)Cr*i  - 

1 1 jJ 


X*  "«  -  V  -  V 

S  mi j^qij(rij)  '  wij-^  ^rij  ’  rij^ 


i  >3 


where  the  second  inequality  follows  from  lemma  A.l. 


OO)  -0(0)  »g:#T”lJ[pU(*lV-Nj]I.*lV  V 

'  0<V*ij 

i^%0P  yniJ  '  Xij-*  ^ij  ‘  Sj^ 

+  y* 

2  Y  nij  LPij^S'j^  "  xip  ^ij  '  S j ^ 

1  J  '  J 


mij  (r*j)  -  trij  -  rij3 


*  rl 

0<r1j<rij 


+  U:r*.=0  Cq<J  tr*J)  '  UlJ]  [r'J  '  rii] 
*  J 


♦  £ 


S  *  d  mi j  ^i  j  ^rij^  "  yij^  ^rij  ”  rij^ 
1,J*r1j  ri j 


Now  from  (3.1)  and  (3.2)  the  first  and  fourth  terms  are  zero.  All  other 
terms  are  greater  than  or  equal  to  zero.  Hence 


0(1)  -  0(0)  =  op  -  Jp*  >0 

That  is,  0  at  any  feasible  point  P  is  greater  than  or  equal  to  0  at  ?* 
This  completes  the  proof  of  sufficiency  of  (3.1)  through  (3.4). 
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APPENDIX  B 


In  this  appendix  we  give  the  proof  of  theorem  3.2,  which  deals  with 
the  convergence  properties  of  the  joint  flow  control  and  routing 
algorithm  H.  Recall  that  the  objective  function  0  is  given  by 

J  ■  E«ik  <F?k>  ♦  S  Gik  (F1k)  *  2>a  Va  Ua)  *  £">a  Ea(ra) 
i,k  i,k  a  a 

Define  the  vector  u  by  the  tuple  u  =  (i,  r,  Fv,  Fd) ,  where  l  is 

the  vector  of  voice  packet  lengths,  r  is  the  vector  of  data  input  rates, 

Fv  is  the  vector  of  link  voice  traffics,  and  Fd  is  the  vector  of  link 

data  traffics.  Let  un  denote  the  vector  u  at  iteration  n.  Assume 

the  algorithm  is  started  with  a  feasible  u°  such  that  J(u°)  <  JQ  , 

for  some  positive  J„.  Let  the  set  U  be  defined  as  follows: 
o  o 


U0  -  {u|J(u)<J0} 


(B.l) 


Then  u  e  U  .  Define  A (9)  as  follows: 


A(9)  =  J(u9)  -  J{u") 


(B.2) 


where  u9  *  8un+^  +  (1  -9)un  ,  0  <  9  <  1.  Then  for  all 


u9  e  UQ,  a(9)  is  continuous  in  9,  and 


A(0)  =  0 

A ( 1  )  =  J(un+1)  -  J(un) 


(B.3) 

(B.4) 


The  Taylor  expansion  of  A{0)  about  zero  is 

Me) 

•  •£  ^  c-r1 

i  1 

n  i 
-  u.  ] 

♦  1  °2 

i  j  1  J 

t-r1  -  “p 

e* 

u=u  : 

0  <  0*  <  0 

Thus 

dA(0) 

d0 

s 

£  Mini  tuyi  .  u» 

] 

0 

=0 

1  1 

3 

s  C8jk  <f™)  *  =;k 

-  Fvn  1 
hik  J 

i.k 

+ 

T'' G*  (Fn  )  rF<^n+^ ) 
^bik  lhikJ  uik 

-F?P 

i  »k 

+ 

2  »;  (*I>  cc1  - 

9 

+ 

E;  (rj)  [r^  - 

d 

Lemma  B.l 

Assume  that  uneUQ.  Then 

^  o. 

0=0 

Proof  From 

(3.25) 

and  (3.26) 

x. 

vn  vn  r.vn 

!  Ap  *  L**  Cdp 

-  dvn(a)]  [S™  - 

$v(n+1 ) j 

Pepa 

*Pv 

(since  e^* 

=  0) 

-86- 


That  is. 


vn  .vn 

e»  l‘r> 


SC  CSp"  -  svp<"+,>]  -  dvn(a)  Sts'"- s'("+1>: 


*  S>  [Svn  -  sv(n+1)]  -  dv?  [s™ 

** up  L  p  p  J  P  D* 

P  V 


-  s^1*] 

P„ 


dVn(a) 


S.  tc-  C"^ 

p^p  p  K 


Now  dv* 
P 

v 


*  dvn(a) ,  and  Sv2  =  Svn(a) 


-  S 


vn  .vn  _ 
e„  a  - 
P  P 


-Sc  -  C  3 


But  d 


*  dvn(a)  CSv(n+1>(a)  -  Svn(a) ] 
*  Ssn(i,k)  [B.k  (F™)  +  Gik  (P^)] 


Therefore, 

S  £  c  c  ■  -£c4<4>  +  4<4>£: 

a  pePa  i,k  a  pePa 

♦  Sdv"U)  tsv(n+1)U)  -  S™U)] 


"St4  <C)  *4  (F?k)]  [F^+1)  -  F™] 


6Y  E"„  dvn(a)  U^1  -  # 


1-hs-vJt 


-87- 


where  the  last  equality  follows  from  (3.15)  and  the  definition  of  Svn(a). 
From  (3.21) 


dvn(a)  =  z;  +  pa(ij) 


"»LE  V  ?  -  -  S  (F-)  ♦  s:k  <F?k>]  F«3 


a  peP. 


i  ,k 


+  By  2>,  h  <0  Cl"+1  ' 


-  «  zj  UJ  -  C’] 


■  -  £  CBik(F-)  +  Gn  (Fjk):  CF'<"+’)  -  F*"] 


i,k 


ik- 


a  na »;  <*!>  t*r'  -  # 


-  6  vSna  z£  [lan  -  C’  3 


(B.5) 


where  the  last  equality  follows  from  the  definition  of  PaUa);  see  (3.9) 
Now  from  (3.22),  |aJ  -  tjj+1 1  <_  nt|z!J|.  Moreover,  when  zj  >  0, 
a"  -  >  0;  and  when  z"  <0,  a"  -  a"+1  <  0.  Thus 

a  a  —  a  —  a  a  — 


i>:  -  o  -  i^i  t‘an  - 


Therefore,  Sy£  n,  zj  [*a  -  i^1]  >  BY  rj1  IX^  -ij+1]2 


and  -  BY^n,  zj  [tj  -»;]<-  BY  n'1  £  -  ff1]2  » 


,n  »n+l t2 
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Thus  from  (B.5) 


a  peP,  i,k 


(•„'  /..vrK  4.  r '  \l  TFv^n+^  -  Fvnl 

CBik  {?i*)  +  Gik  (Fik)]  [Fik  FikJ 


-  Z 


V  na  »;  (!»)  [  r1  -  *!l 


-  8  Y  n'1  Z"a  CC'  -  1"J2 


(B.6) 


Also  from  (3.26)  nv  e™  .  Thus  A^n  e™  >_  r^1  [A™] 


Therefore,  from  (B.6)  we  have  that 


«;'E  Z  i>p"]2  i  -  Z  c»;k  cf«)  *  s:k  (p;k):  cf* 

a  peP.  i >  k 


v(n+l) 

ik 


Fvnl 

hikJ 


*  "a  \  t*S»  K'  -  O 


-  6  '  n'1  Z  "a  Cr  -  *a"]2 


(B.6a) 


Or 


T,  cb! 


i  ,k 


(FlC)  +  (F"  )] 


ik  ''ik'  ik  '  ik' J  L'ik 


[Fr'n+1)  -  fiu.S 


If  •  Tvna  »;  ttJJ) 

a 


<  -  &  Yn 


''In,  [if1 


i  a 


a  '■  a 


-  O2  -  c'E  EcO2 


a  peP, 


Similarly,  it  can  be  shown  that 


£  ^  (F*,)  [F^n+1  >  -  fJ“1  ♦  £  *a  £>?)  [rf1  -  ra"  ] 


i  »k 


ikJ  "a  “a'1  a'  ua 

a 


i  -  "‘J  S  ma  C^1  -  rna]  -  n'1  £  £  [*f  ]2 

<& 


Therefore, 


dA(6) 


de 


1  -  n-'sir£  na  [C1  -  2 


0=  0 


-  %'£  £  tO2  -  "j'£  £  [if]2 


v  »  p3Sa  o 


a  peP. 


-  n 


£  ma  tC  -  2 


(B.7) 


<  0 

This  completes  the  proof  of  lemma  B.l. 
Lemma  B.2.  Assume  un  ell0.  Then  ^9) 


<  0 


8=0 


if  u  is  not  an  optimal  point. 


Proof  From  the  result  of  lemma  B.l  we  have  that  ^9) 


<  0.  To 


9=0 


show  that  the  inequality  is  strict  it  suffices  to  show  that  if  un  is  not 


an  optimal  point  then  at  least  one  of  the  sums  in  (B.7)  is  strictly 
positive.  Assume  that  none  of  the  four  sums  in  (B.7)  is  strictly  positive. 
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Then  for  voice  we  have  the  following  situation: 

(i)  _  2.”  .  If  0<  £? <  a. .  ,  then  this  situation  implies  that 

a  -  d  a  IJ 

Z "  3  dvn (a)  -  o  (*")  =  0,  and  the  marginal  gain  in  voice  quality 

d  a  o 

is  equal  to  the  marginal  voice  cost  of  congestion.  If  l"!  -  0, 

a 

then  z"  >0  and  dvn(a)  >  paUj).  And  if  l1}  3  a..,  then  Z?<0 
a  a  a  a  ij  a 

and  dan<Pa(^). 

(ii)  Avn  3  0,  which  means  that  e^n  3  0  since  by  assumption  S^n > 0 . 

P  P  P 

That  is,  e^11  3  d^  -  dvr\a)=  0,  and  all  voice  traffic  flows  on 
shortest  paths. 

By  a  similar  argument  for  data  we  have  that 

!=  ddn(a)  if  0  <  rj  <  rd 

>ddn(a)  ^  ra  *  ra 

<  ddn{a)  if  r!J  3  0 

(iv)  All  data  traffic  flows  on  shortest  paths. 

But  results  (i)  through  (iv)  are  precisely  the  Kuhn-Tucker  conditions 
for  optimality  given  in  (3.1)  through  (3.4).  Thus  by  the  sufficiency 
part  of  theorem  3.1,  proved  in  Appendix  A,  un  is  an  optimal  point. 
Therefore,  if  un  is  not  an  optimal  point  at  least  one  sum  in  (B.7)  must 
be  strictly  positive,  thus  proving  the  lemma. 

The  implication  of  lemma  B.2  is  that  for  small  enough  9,  A (9)<  0. 
Since  A(0)  is  continuous  in  0  for  u  eU0,  either  A(9)  <  0  for  all 


I 


6  <  1  or  there  exists  a  1  such  that  &(e)  <  0  for  e  <  8-j  and 
a(0-|)  *  0.  That  is,  the  plot  of  A(e)  takes  the  form  of  either  curve  (i) 
or  curve  (ii)  of  Fig.  B.l.  Let  be  the  smallest  e,  0  <  3  <_1 ,  for 

which  a (e )  =  0.  If  no  such  e  exists,  then  let  a (0, )  =  A ( 1 )  • 


Figure  B.l  Possible  Forms  of  A (9) 


Assume  that  there  exists  e-|  <  "* »  such  thatA  ($1  )  =  0.  The  Taylor 
expansion  of  A  (8,)  (about  zero)  is 
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From  lemma  B.l 


3J(un)  r,n+! 
3ui  LUi 


dA(9) 

de 


9 


=  0 


9  <  9] 


0 
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The  last  equality  follows  from  the  fact  that  we  are  taking  partial  and 
not  total  derivatives.  By  the  same  argument  we  have  that  all  terms 


2  2 
fiu^uT'  *  which  ui  uj*  are  zero  except  that  t  0-  Thus 

i  j  3Fik3Fik 


^2  [U?+1  -  u?]  [u^1  -  ulll 

j  3ui  3uj  L  1  L  j  jj 


a' 

U  =  U° 


(F?k>  +  l4  (Fik  ”  C^"+1)  -  O2 

i,k 


*  E  <k  (<)  c^r’  -  ^2 


-  2  E  «;»  (F?L»  tF^+»  -  F™]  [F?<"+’>  - 


ikJ  L  ik 


+  E  T  "a  VI  <*1>  CC1  • 

a  * 


+  E  "a  <  (f f )  Cf"+1  -  rtf 


(B.8) 


SinceA  (0)  is  continuous  in  0  for  u9eU0,  the  second  derivatives  exist 
over  the  compact  region  0  <_  0  <_  9^. 


Lemma  B.3  Let  M  *  max  max 


u:  J(u)<JQ 


MO 


for  some  positive  0Q.  Then 


E"avl  <4’>  [C’-o2 1  \  E  %  ur1 


Proof  Recall  that  UQ  =  (ujJ(u)  <_ JQ}  for  some  positive  JQ.  Since 
we  assume  u6  e  U.  for  9c[0,l],  then  V_( A®  )  <  M0.  The  proof  of  the 

u  a  a  —  X» 

lemma  follows  then  from  the  assumption  of  the  lemma. 


Lemma  B.4.  Let  M_  =  max  max 
-  r 


Then 


E"  (rj  for  some  positive  J  . 
a  u:J(u)<Jo  a  a  0 


Proof  Similar  to  the  proof  of  lemma  B.3. 


Lenwa_ B.5  Let  M  *  max  max  [b"  (fY.  )  +  G"  (F..)] 

v  i,k  u:J(u)<J  1lc  lk  1k 

for  some  JQ.  Furthermore,  let  the  cardinalities  of  the  sets  Ay  and  Pfl 

be  1 Avt  and  |Pj  ,  respectively.  Then 


£  (fjh')  *  o;k  (<>]  -  o2 


<  Hv  N2|A  l2|p  |3  £  £  [a”n]2 

v  v  a  a  peP  p 


Proof  For  the  same  reason  given  in  the  proof  of  lemma  B.3,  we  have 

I  I 

‘hat  B^F*  )  *  <k  (F»k)  ‘  Hy.  Thus 

5  :^n+,)  -  Q2i\  Ek  cf^’>-  02 


From  (3.15)  we  have 

rcv(n+!)  _  pVn-]2  _ 
uik  MkJ 


that 

[E  E  «p  (i.w  [s^+’>  - 

a  PcPg  k  h  H 

I  \l  ]Pai  E  E  5p(1.k)  [s“(n+,)  -  S™]2 

d 


< 
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where  the  inequality  is  the  Cauchy-Schwarz  inequality.  Since  we  make 
incremental  changes  at  each  iteration,  the  largest  change  in  any  path 
flow  cannot  exceed  the  sum  of  all  changes  in  path  flows  for  the  same 
source  -  destination  pair.  Therefore, 


ls;(n+,)  -  s£nl  I  ™  .  p,pa 

ThUS  [Sj(n+1)  -  S^32  <  C22  a*"]2 

pePa 

<  iPa!  T,  [Ap"12 

pePa 

where  the  last  inequality  is  the  Cauchy-Schwarz  inequality.  Therefore, 


<  |AV12  |Pa|3  £  £  5P(1'k)  ^ 


a  peP. 


That  is,  [F^n+1)  -  F^]2  <  |Ay|2  |Pj3  23  23  M  MO  [A™]2 


a  peP. 


And 


E  [f?<"+1)  -  f™:2<n2  uvi2  iPai3  E  E  tAp"]2 


i  ,k 


a  peP. 


Multiplying  both  sides  by  My  completes  the  proof  of  the  lenma. 


Lemma  B.6  Let  M.  *  max  max  G..  (F.  J 

d  i,k  u: J(u)<J  1k  1k 
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for  some  J_.  Furthermore  let  |  Pa  |  be  as  defined  in  lemma  B.5,  and  let 

0  a 


[ Ad|  be  the  cardinality  of  the  set  A^.  Then 


E  4  (F^  [4"+,)  - 4]2  n2  m2  iM3  E  E  cO: 

i  ,k  a  PeP, 


Proof  The  proof  is  similar  to  the  proof  of  lemma  B.5. 


Lemma  B.7  Let  and  |Ad|  be  as  defined  in  lemma  B.6,  and  let 


[ A  | and |  P  I  be  as  def i ned  i n  1 emma  B.5. 

V  a 


Then 


i,k 


4  <  Cffi"*1 }  -  F*] 


ikJ 


i  ?"dN2  l\l2lpal3  E  E 

a  pePa 


Md  n2I  Ad  1 2 1  pal  3  E  E  Cip"]: 

a  pcP. 


Proof  We  know  that  (x^  -  y^)2  =  (x2  +  y2  -  2x..  y.)  >_  0 

i  i 

£xiyi  <7  +?2>?  (B-9) 


Thus 


i  1  i 

Now  from  the  assumptions  of  the  lemma, 

E  4  <4>  [4"+,)  -  to  ^ 


i,k 


1  Hh  E  <nt'>  ■  TO  -  Ffr 


'*  t$  ”'k 


ik- 


P 
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If  we  let  x.j  =  F^n+1  ^ 


cvn 

hik 


and 


=  F<?(n+1) 


i  'ik 


-  F( 


dn 

ik 


then  from  (B.9)  we  obtain 


52  4  <4>  r4n+1)  -  O  c4n+,)  -  43 


i  ,k 


<  i  M  T  rFv(n+1)  -  FVnl2 

±  7  "d  f-£  Lhik  hikJ 


*?"d  -F*]2 

1  ,  K 

"d"2  IM2  lpal3  52  52  [a™]2 

a  PePa  p 

+  JMdN2  lAdl2  lpal3  52  23  [Ap"]2 

where  the  last  inequality  follows  from  lemmas  B.5  and  B.6. 

From  the  results  of  lemma  B.l,  and  lemmas  B.3  through  B.7  we 
obtain 

A(8,)  18,^8, (Mv  Md)  N2  |AV|2  |Pa|3  -  rf,1]  r  52  k? 

a  pePa  K 


*  8^8,  Md  N  [Ad|‘  |Pa|J  - 


-  52 


dn.2 


] 


a  PePa 
.n+1  „n,2 


*  V4BY]En>U8+1  -i"a] 

a 

+  9t  [7  0!  Mr  -  n;1]  5Zma  [rj+1  +  r£]2 


Lemma  B.8  Let  the  scale  factors  be  chosen  as  follows: 


nv  *  tj  (Mv  +  Md)  N2  |AV|2  [Pjl3]'1 
nd  *  CMd  N2  |a/  |Pa|3r' 

n*  s  s  V 

Ar  *  [j  Hr]”' 

Then  if  un  is  not  an  optimal  point,  A(6^)  <  0. 

Proof  Substituting  the  n's,  as  defined  in  the  lemma,  in  (B. 10)  we 
obtain 

^(6^)  <  01(MV  +  Md)  N2  |Ayi2  IPJ3  [•,-!]£  X)  CApn]2 

a  pePa 

+  ^  Md  N2  !Ad|2  |  P  i3  [9]  -  1]2  S  Op”]2 
1  d  a  3  1  a  peP  p 

a 

+  1  9!  C91  -  1]  £  "a  CC1  - 

a 

+  r  8,  »r  ce,  -VjEst'T1  -^2 

a 

Since  we  assume  9^  <_  1 ,  we  have  that  A (8 ^ )  <_0.To  show  that  A(9-|)  <  0 

I 

if  the  n  s  are  selected  as  stated  in  the  lemma,  it  suffices  to  show 
that  at  least  one  of  the  sums  above  is  strictly  nonzero  if  un  is  not 
an  optimal  point.  But  this  is  what  we  proved  in  lemma  B.l.  Thus  the 
proof  of  the  lemma  is  completed. 


* 
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The  result  of  this  lemma  shows  that  there  exists  a  choice  of 
scale  factors  such  that  A(0^)<  0.  Therefore,  with  the  scale  factors 
chosen  as  in  lemma  B.8,  the  assumption  that  A(e -j )  =  0  for  0]  <  1 
is  contradicted.  Thus  if  un  is  not  an  optimal  point  and  the  scale 
factors  are  chosen  as  in  lemma  B.8,  a(0-j )  <  0  for  0  <  0-j  <_  1. 

And  in  particular,  a(1)  <  0.  That  is,  J(nn+1)  -  J(un)  <  0.  By 
induction,  J(un)  <  J(nn  ^)  <  .  .  .  <  J(u°).  But  we  started  with  a 
feasible  u°  such  that  J(u°)  <  JQ.  Therefore  J(un)  <  JQ.  Thus  the 
algorithm  is  a  descent  algorithm;  and  if  u°  eUQ,  the  u11  eU0  for 
all  n  >_  1 . 

Lemma  B.9  The  mapping  H  is  closed. 

Proof  We  define  closedness  of  a  point- to-set  mapping  in  the  same 
way  it  is  defined  by  Luenberger  [25].  That  is,  let  H  be  a  point- to- 
set  mapping  which  maps  points  in  a  space  X  into  subsets  of  another 
space  Y.  Then  H  is  said  to  be  closed  at  x  cX  if  the  assumptions  that 

(i)  {xn>  -►  x  ,  for  xn  e  X 

(ii)  {yn}  -*•  y  ,  for  yn  e  H(xn) 

imply  that  y  eH(x). 

Let  {wn}  *  {Spn  ,  sjjn  ,  ,  rj}  be  a  sequence 

which  converges  to  w  =  (Sp,  sjj,  $,a,  ra).  Then 
s»n(a)  .  £  $vn  ,sv(j) 
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Si  nee  fJJ 


and  since  S^n  -*■  F^  if  we  started  with 

u°  eUQ  (since  by  the  previous  result  all  subsequent  un  are  also 


vn 

P 

vn 


elements  of  UQ). 


Thus  d 


vn 


2  ^4  <Fik> +  <4  “ 

(i.k)ep 


B.,  and  G..  are  continuous  in  fV!\  And 


ik 


ik 


ik' 


since 


dvn(a)  =  Min  dY”  ^ dv(a) 

Pepa  » 


.vn  .  dvn  .  dvn(a) 


'P  -P  -  CP  * 

Therefore,  =  min  [sjn  ,  e™]-*  a* 


Let  the  sequence  {  w  }  =  {  sVl  $j!n  ,  i",  r"}  be  another 

n  p  p  a  a 

sequence  that  converges  to  w  =  (§)!,  S^,  £  ,  rj .  Furthermore,  let 

p  p  a  a 

the  sequences  (wn)  and  (w  }  be  related  as  follows: 


Sp  -  Ap 


if  p  t  P!n 


:vn 


Svn(i) 


-§*n 


P^P 


,-vn  ,vnx  *n 

(Sp  -  Ap  )  if  P  »  Pv 


where  py  is  an  optimal  voice  path  at  iteration  n. 


'dn 


edn  „  dn 
sp  •  4p 


Sdn(a)  - 


If  p  t  pd" 


p  ,  -  W  "P-P* 


-101- 


Jrn 

where  pd  is  an  optimal  data  path  at  iteration  n. 


7" 

n*Za 


if  0  <  *  ^  Za  <  aij 


‘U 


if  An  - 
d 


\  zn  >  a. . 
*£  a  - 


if  A 


n  Zn  <  0 
A  a  — 


.  | 

‘  l 

where  z”  and  x"  are  as  previously  defined.  Note  that  Z1?  and  x"  converge 

da  a  a 

to  Z,  and  x,,  respectively,  since  these  terms  (i.e.  Z ”  and  x")  are 

d  a  a  a 

continuous  in  wn. 

Now  we  have  that  Spn  -  Apn  converges  to  Sp  -  Ap  for  all  p  and, 

^y  Aw  y  y 

by  assumption,  Sp  converges  to  Sp  for  all  p.  Suppose  Sp  t  Sp  -  Ap 

for  some  p.  Then  there  exists  some  e  >  0  such  that  |S  -  (Sp  -  Ap)]=  £• 

★n 

Also  there  exists  some  nQ  such  that  for  n  >_  nQ  and  p  f  pv  , 
lspn-  (Spn  -  Ap") |  <  e/2.  Thus  for  n  >  nQ,  p  -  p*n  .  That  is. 


ra  ~  ^  x!)  if  0  <  r"  -  r\  x,  <  r^ 
did  a  r  d  d 


..  n  n  .  d 

’f  r  -n„  x  >  r 
a  r  a  —  a 


if  rn  -  n  xn  <  0 


Svn  +  sv 
P  P 


Sp  -  *p  "P*pJ 


Sv  (a) 


y:  (Sp  -  Ap)  if  ?  * 

p*p! 
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where  py  is  some  limiting  optimal  path.  By  the  definition  of  closed¬ 
ness  [25],  we  conclude  that  H  is  closed  at  Sp  since  Sp  eH(Sp). 

By  a  similar  argument  we  also  conclude  that  H  is  closed  at  S^,  since 

A  J  J 

eH(Sp).  Since  Z"  converges  to  -n^  Zfl  for  all  a,  and 

Ah  A 

since  by  assumption  i,  converges  to  £  .  we  can  apply  the  above 

q  a 

argument  to  conclude  that 


^a  "  Za 


2,  -*•  z  -  J  a .  - 

a  a  \ 


1f  0  <  “a  -  "t  2a  ‘  a1j 

1f  la  '  Va  i  aij 


ff  *a  '  Va  i  0 


Thus  H  is  closed  at  £  since  £.  eH(£  ).  And  by  a  similar  argument 

d  a  a 

A 

we  conclude  that  H  is  also  closed  at  r_,  since  r  £H(r  ).  Thus 

a  a  a 

A 

H  is  closed  at  w  since  w  eH(w),  and  this  completes  the  proof  of 
the  lemma. 


Proof  of  Theorem  3.2  Let  U  be  the  following  compact  set: 

(J  =  (ujJ(u)  u  is  a  feasible  point  of  (2.14)  through  (2.16)} 

Then  H  is  a  closed  point-to-set  mapping  of  points  in  U  into  sets  of 
points  in  U.  Also  by  lemma  B.8,  the  sequence  { J (un) }  is  a  strictly 
decreasing  sequence.  That  is,  if  un  is  not  an  optimal  point  then 
J(n°  )  <  J(un),  where  un+^  eH(un).  Therefore,  by  the  global  convergence 
theorem  [25,  p.  125],  the  limit  of  any  convergent  subsequence  of  the 
sequence  {un}  generated  by  H  is  a  solution  to  (2.13).  This  completes 


the  proof  of  theorem  3.2. 


Proof  of  Necessity  of  Theorem  3.1 


From  lemma  B.8,  we  observe  that  0  strictly  decreases  if 
equations  (3.1)  through  (3.4)  are  not  satisfied.  Thus  at  any  point  u 
that  these  equations  are  not  satisfied  the  algorithm  produces  a  better 
point,  which  means  that  u  cannot  be  an  optimal  point.  This  then 
proves  the  necessity  part  of  theorem  3.1. 
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